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1.  Research  Objective* 

For  the  period  Dec.  15,  1983  -  Dec.  14,  1984,  the  proposed  work 
plan  included  three  tasks: 

(a)  To  complete,  develop  and  verify  a  quasi-two  dimensional  MPD 
channel  model  which  was  under  development  at  proposal  time.  This  model 
included  many  physical  effects  which  had  not  been  previously  addressed 
(finite  ionization  kinetics,  particle  interdiffusion,  etc.),  and  was 
specifically  constructed  in  an  attempt  to  verify  and  refine  the  onset 
predictions  of  Ref.  1. 

(b)  To  conceive  and  analyze  several  design  variations  including 
axial  fields,  segmented  electrode  construction  and  contoured  wall 
designs.  This  was  to  be  done  by  and  large  analytically,  with 
initiation  of  a  numerical  model  development  also  contemplated. 

(c)  Definition  and  design  of  a  small-scale  experiment  to  test  one 
of  the  above  concepts,  to  be  performed  on  an  outside  facility  with 
support  from  MIT. 

(d)  To  coaq>lete  an  ongoing  study  of  the  local  stability  of  an  MPD 
plasma,  including  both  the  partially  ionized  and  the  fully  ionized 
range  of  operation,  and  assessing  the  effects  of  finite  size,  heat 
conduction,  species  diffusion  and  finite  self-magnetic  field  effects. 

The  following  section  gives  an  overview  of  the  work  accomplished 
in  each  of  these  tasks,  and  aakes  reference  to  several  Appendices, 


where  more  detailed  results  and  discussions  are  contained. 


The  bulk  of  the  work  proposed  under  Task  (a)  above  was  in  fact 
completed  during  the  first  year  of  the  Grant,  and  reported  in  our  1st 
Yearly  Report.  The  work  showed  the  effects  of  geometrical  parameters 
on  degree  of  ionization,  wall  recombination  losses,  etc.  One  very 
important  result  was  the  fact  that  partial  ionization  had  only  a  minor 
impact  on  the  predicted  "onset"  current.  We  have  used  this  result  to 
justify  the  continued  use  of  simplified  full-ionization  plasma  models 
designed  to  illuminate  other  effects,  such  as  those  of  a  predominantly 
dynamic  nature.  It  is  clear,  however,  that  any  assessment  of  losses  or 
efficiency  must  be  based  on  more  complete  plasma  formulations. 

The  original  Baksht  model  (Ref.  1),  as  well  as  our  own  extension 
discussed  above,  were  specifically  confined  to  the  "onset”  condition. 
This  was  done  by  imposing  a  condition  of  vanishing  plasma  density  at 
the  anode  wall.  It  is  of  interest  to  inquire  about  the  approach  to 
onset  from  below,  since  that  may  provide  guidance  as  to  how  close  to 
that  limiting  condition  it  is  safe  to  operate,  and,  in  addition,  it 
should  illuminate  the  physical  events  responsible  for  or  triggered  by 
the  gradual  depletion  of  the  anode,  such  as  the  growth  of  the  anode 
voltage  drop.  We  report  such  an  analysis  in  Appendix  1  to  this  Report. 
The  conclusions  are  that: 

(a)  Onset  symptoms  are  absent  up  to  currents  90k  of  critical 

(b)  Anode  drops  remain  close  to  those  corresponding  to  the 
floating  potential  for  the  whole  channel  also  up  to  about  90k  of  I" . 

(c)  Anode  plasma  density  decreases  with  current  at  about  the  same 


A*  v!‘. 
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rata  as  cathode  density  up  to  near  P ,  then  drops  rapidly  to  near  zero. 
The  anode  drop  is  negative,  except  very  near  onset. 

In  work  related  to  Task  (b)  above,  we  undertook  local  analyses  of 
various  parts  of  a  typical  MPD  thruster,  in  an  effort  to  isolate  and 
illuainate  the  dominant  physical  effects  in  each. 

Near  the  inlet,  where  the  flow  speed  is  still  low,  the  Bain 
component  of  the  local  electric  field  is  purely  ohaic,  producing  the 
strong  local  dissipation  which  leads  to  rapid  ionization.  One  problea 
which  has  recurred  in  eost  channel  designs  is  a  strong  current 
concentration  at  the  cathode  root,  with  erosion  of  both  the  cathode  and 
the  backplate.  While  one  part  of  this  has  been  shown  in  the  previous 
analysis  (App.  1)  to  be  associated  with  the  reduced  uB  back-oaf  in  that 
region  (and  becoaes  especially  strong  at  onset ) ,  we  expect  one 
additional  component  of  current  concentration  at  the  cathode  root  due 
to  the  peculiarly  skewed  current  pattern  occurring  in  the  presence  of  a 
finite  Hall  effect.  To  analyze  this  effect,  we  start  in  Appendix  2  by 
formulating  a  general  tensor  Oha’s  law  in  a  fora  which  shows  a  clear 
division  of  the  local  field  into  the  gradient  of  a  generalized 
potential,  plus  a  generalized  back  eaf  (these  generalizations  arise 
froa  the  Hall  and  electron  diffusion  effects).  This  leads  naturally  to 
an  induction  equation,  governing  the  magnetic  field.  We  then 
particularize  to  saall  flow  velocities,  and  show  that  the  field  B  obeys 
v*B  =  0,  despite  the  Hall  paraaeter  gradients  oc curing  due  to  density 
gradients  (which  are  theaselves  due  to  axial  currents).  It  is  also 
shown  that  in  this  region,  and  for  a  fully  ionized  plasaa,  electron 


diffusion  has  the  effect  of  reducing  by  a  factor  of  2  the  effective 
value  of  the  local  Hall  parameter.  By  Baking  the  additional 
sieplification  of  taking  a  Bean,  constant  Hall  pareBeter,  we  are  able 
to  obtain  analytical  solutions  for  the  current  density  distribution  on 
the  electrodes,  including  a  singularity  of  the  inverse  factional-power 
type  at  the  cathode  root  (and  a  siailar  zero  at  the  anode  root). 

The  identification  of  anode  depletion,  driven  by  axial  currents, 
as  the  aost  likely  trigger  for  onset,  had  led  us  to  propose  (Ref.  3)  a 
segaented  channel  construction,  with  only  two  external  terainals  due  to 
diagonal  cross-connection,  designed  such  as  to  cancel  the  axial 
currents  by  the  establishaent  of  appropriate  axial  electric  fields. 

There  are  aany  questions  to  be  answered  regarding  this  design, 
which,  however,  holds  the  potential  for  significantly  delaying  onset, 
and  hence  increasing  thruster  power  and  efficiency.  In  Appendix  3  we 
have  provided  a  siaplified  analysis  to  show  how,  under  conditions  where 
axial  current  has  indeed  been  suppressed,  the  transverse  current  and 
the  voltage  are  distributed  in  the  axial  direction.  This  leads  to  the 
establishaent  of  liaitations  for  this  type  of  thruster:  essentially, 
the  ratio  of  current  to  aass  flow  aust  be  above  a  certain  liait  (shown 
to  correspond  to  a  Hall  paraaeter  of  0.6  for  Argon),  or  else  not  enough 
voltage  overlap  develops  between  the  anode  and  cathode  walls  to  allow 
cross-connection.  It  is  shown  that  thrusters  with  a  less  than  about  2 
g/sec  can  easily  be  designed  to  avoid  this  liaitation,  while  still 
reaaining  below  traditional  onset  liaits  (and,  in  addition,  the  new 
onset  liait  is  now  likely  to  be  higher). 
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;  Work  has  also  baan  initiated  on  a  nodal  of  the  effect  of  wall 

i 

contouring  on  the  diatribution  of  current  and  other  quantities.  This 
model  exploits  the  simplifications  afforded  by  the  high  values  of  the 
typical  magnetic  Reynolds  number  based  on  length  (5-10),  and  of  the 
1  ratio  of  magnetic  to  gas  pressure.  Results  are  expected  within  the 

i 

next  reporting  period. 

The  experimental  work  has  been  deferred  pending  acquisition  of 
advanced  instrumentation  to  be  acquired  under  a  parallel  DoD 
Instrumentation  Grant  and  completion  of  our  contoured  wall  analysis. 
It  is  our  intention  to  design  a  contoured  channel  for  even  current 
distribution  (minimizing  hot  spots)  and  to  conduct  detailed  diagnostic 
tests  comparing  its  performance  and  flow  and  current  patterns  to  theory 
and  to  conventional  thrusters.  This  will  also  serve  to  assess  the 
feasibility  of  a  later  segmented  channel,  a a  discussed  above.  Tests 
will  be  at  an  outside  facility  to  be  selected,  using  an  MIT  provided 
thruster  and  also  MIT  optical  diagnostics  and  digital  data  acquisition 
equipment. 

Finally,  Appendix  4  gives  a  complete  account  of  our  investigations 
on  the  local  stability  of  an  MPD  plasma,  a  part  of  which  was  advanced 
in  our  previous  Report  (Ref.  2).  The  new  work  consisted  of  a 
clarification  of  the  types  of  potential  instability  modes  depending  on 
the  operating  regime,  mainly  the  degree  of  ionization  a.  It  was  found 
that,  up  to  a  -  0.9  the  predominant  effect  is  the  classical  ionization 
instability,  involving  mainly  variations  in  electron  density,  whereas 
very  near  a  =  1,  a  purely  static  node  involving  mainly  tesperature 
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variations  becomes  predominant.  Heat  conduction  stabilizes  this  node 
for  Moderate  Hall  parameters,  but  this  effect  is  weakened  and 
eventually  disappears  at  a  critical  value  of  this  Hall  parameter.  This 
is  qualitatively  similar,  but  for  a  completely  different  reason,  to  the 
role  played  by  the  Hall  parameter  in  catalyzing  the  ionization 
instability.  For  MPD  conditions,  it  is  concluded  that  striations  with 
wavelenghts  below  some  3  cm  will  be  well  damped;  since  regions  of 
near-complete  ionization  are  unlikely  to  be  much  more  than  this  size, 
the  expectation  is  that  the  static  node  will  not  became  unstable. 
Again,  this  conclusion  is  similar  to  that  reached  earlier  (Ref.  2)  on 
the  sense  that  the  ionization  mode  is  unlikely  to  be  unstable  due  to 
the  small  special  extent  of  high  Hall  parameter  regions.  Taken 
together,  these  results  tend  to  discound  the  possibility  of  the  onset 
phenomena  being  directly  associated  with  bulk  plasma  instability,  and 
to  reinforce  the  anode  depletion  viewpoint. 

3.  Personnel  Associated  with  the  Research  Effort 

The  research  was  led  by  Associate  Prof.  Manuel  Mart inez-Sanchez , 
with  the  cooperation  of  graduate  Research  Assistants  D.  Heimerdinger 
(Doctoral  Candidate)  and  Tze-Wing  Poon  (Masters  Candidate) .  A  new 
Research  Assistant,  J.  Marc  Chanty  joined  the  group  towards  the  end  of 
the  period,  and  will  be  working  on  numerical  MPD  simulation. 

1.  Baksht,  F.Q.,  Moizhes,  B.Ya  and  Rybakov,  A.V.,  "Critical  Regime 
of  Plasma  Accelerator”,  Zhurnal  Tekhnicheskoi  Fiziki,  Vol.  43, 
pp.  2568-2573,  Dec.  1973. 
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Quasi  two-diaensional  MPD  flow  theory  -  A  generalization  of  the  theory 
of  Baksht.  Moizhea  and  Rybakov  to  conditions  below  onset 

Introduction 

The  siaple  theoretical  treataent  of  MPD  flow  given  in  Ref.  1  is 
reaarkably  successful  in  predicting  the  basic  trends  and  even  the 
correct  Magnitudes  of  conditions  at  the  onset  of  instability  and 
erosion.  Additionally,  this  provides  a  clear  Mechanistic  description 
of  the  reasons  for  the  observed  "onset"  phenoaena,  to  be  contrasted 
with  the  soaewhat  obscure  arguaents  contained  in  alternative  theories, 
such  as  those  based  on  the  "critical  ionization  velocity"  or  the 
"ainiaua  voltage  principle". 

In  this  section  we  will  present  a  theory  that  essentially  accepts 
all  of  the  siaplifications  introduced  in  Ref.  1  for  the  treataent  of 
the  flow,  but  reaoves  the  restriction  as  to  the  anode  boundary 
conditions,  so  that  the  results  can  span  a  range  of  currents  below  that 
for  which  anode  depletion  occurs.  This  will  provide  additional  insight 
on  at  least  three  aspects  of  the  problea:  (a)  the  variation  of  anode 
voltage  drop  with  current,  (b)  the  variation  of  overall  transverse 
voltage  with  current,  and  (c)  the  progressive  decrease  of  near-anode 
density  as  the  critical  (onset)  condition  is  approached. 


2.  Foraulation 


12 


The  basic  assumptions  are  as  follows: 

(a)  Two-dimensional,  rectangular,  slender  channel  »  xr) 

ay  ax. 

(b)  Near  full  ionization  (i^  =  nj  =  n,  n,  <<  n, ) 

(c)  Constant  temperature  and  conductivity 

(d)  Axial  Lorentz  force  dominant  over  axial  pressure  gradient 

(e)  Transverse  Lorentz  force  balanced  by  transverse  pressure  gradient 
(weak  transverse  flows) 

In  addition,  we  will  impose  anode  boundary  conditions  similar  to 
those  applying  to  a  Langmuir  probe  in  contact  with  the  local  plasma 
density.  This  will  relate  anode  voltage  drop  to  local  current  density. 
Notice  at  this  point  that  Ref.  1  used  only  a  limiting  form  of  this 
conditions,  namely,  that  at  onset  the  anode  does  not  retard  electrons 
(zero  or  positive  anode  drop).  We  will  find  that  over  most  other 
operating  conditions  the  anode  drop  is  negative,  i.e. ,  the  anode 
rejects  a  fraction  of  the  random  current  in  order  not  to  exceed  the 
current  demands  of  the  channel. 

The  geometry  is  depicted  in  Fig.  1,  with  B  coming  out  of  the  plane 
of  the  paper. 
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Two  types  of  equations  will  be  used:  (a)  axial  balance  equations  for 
which  transversal ly  averaged  quantities,  (  )  will  be  used,  and  (b) 
transverse  balance  equations,  in  which  axial  derivativea  will  be 
neglected.  The  axial  equations  are 


Continuity: 

Moaentua: 


u  A  = 


u  .  dp  _  t 


Mo  *  dx 


(1) 

(2) 


where  ^  is  the  ion  aass,  a  the  aass  flow  and  A  =  D  x  H  is  the 

cross-sect ional  area.  In  Bq.  (2)  we  will  neglect  the  dp  tern  coapared 

dx 

to  the  j,  B  tern  froa  here  on,  and  we  have  used  Aapere's  law  to  relate 
j,  to  B.  An  integration  of  Bq.  (2)  with  the  approxiaate  initial 
condition  u  =  0  gives 


-  B(o)» 
"  2Mo 


(3) 


where  B(o),  the  initial  aagnetic  field  is  related  to  I,  the  total 
channel  current,  through 


I 


(4) 


The  transverse  aoaentus  equation  is 
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—  j.  R  ;  -  2kT  — — 

* 8  dy  ZkT  dy 

(5) 

and,  since  <<  J—  ,  the  current  conservation  equation  reduces  to 

-  0  ;  j,  -  j,  (x) 

(6) 

The  tensor  Oha’s  Ism,  including  electron  diffusion  is 

o  (B  +  ixB^  JS)  *3  +  3x2 

(7) 

g 

so  here  ft  =  0  -  ,  and  P  is  the  Hall  parameter,  which  varies  with  plasaa 

density  as 

P  =  —  B 
en 

(8) 

Solving  Bq.  (7)  for  the  current  c opponents,  and 

approximate  condition  ^  s  0  (to  be  verified  later),  plus 

imposing  the 

din  n  „  Q 
dx  "  U’ 

1  +  &  .  .  .  kT  din  nt. 

-r  Jl  =  (By  ■  +  e  dy  )fi 

(9) 

1  +  &  .  _  n  ,  kT  din  n 

~a  *  =  Sy  -  u»  *  e  ay 

(10) 

one  of  which  can  be  replaced  by 

obtain, 


which  integrates  to 

n*  =  y  +  const.  (13) 

In  order  to  deteraine  the  constant  in  (13)  we  will  iapose  the 
condition  that  the  anode  current  density  equals  the  difference  between 
the  arriving  (retarded)  electron  current  and  the  Boha-liaited  arriving 
ion  current.  This  presupposes  a  negative  anode  drop,  i.e.,  a 
transverse  voltage  profile  of  the  type  shown  in  Fig.  2,  and  can  be 
written  as: 

edVa 

j,  =  en(o)  [  I®-  e  kT  -  v,  ]  (14) 

where  c,  =  V  8kT,  /(m^  )  and  =  V  kTa /■)  .  For  cases  where  dVa  turns 
out  to  be  positive  (only  very  nesr  onset  conditions),  we  will  replace 
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3*  =  «>(o)Si 


(15) 


and  neglect  the  avail  retarded  ion  current.  Both  (14)  and  (15)  can  be 
written  together  aa 


jy  =  en(o)  y 


(16) 


where  d1  ii  a  function  of  dV,  whose  definition  is  obvious  froa  Bqs. 
(14)  and  (15).  Solving  (16)  for  n(o),  we  can  use  this  as  an  initial 
condition  in  Bq.  (13),  with  the  result 


n(y)  *  V  (~)*  +  y 


(17) 


For  use  in  the  axial  equations,  we  need  the  transversally  averaged 
density.  This  can  be  explicitly  found  froa  (17)  as 


n  =  §  -  {  (  (^/  +  -  (f^-)1 1  }  (18) 


Turning  now  to  the  y-coaponent  of  Oha's  law  (Bq.  10),  we  notice 
first  of  all  that 


d*  -  oB*  -r  2  kT  din  n 
-  J»  '  Jf  =  ~ 


e  dy 


(19) 


where  we  have 


use  of  8qs.  (8)  and  (12).  Bq.  (10)  then  reduces  to 


and,  If  the  cathode  potential  is  taken  to  be  zero,  we  can  integrate  to 
obtain  the  voltage  across  the  plasaa  core  as  (with  reference  to  Fig.  2) 

V  -  dV.  -  dVc  =  (J-  +  Z  B)  d  jp  In  (21) 

or 

V  =  dV.  +  dVc  +  (^-  +  u  B)  +  |e  ln(1  *  laf  ~--C~  d)  (22) 

Ji 

Bqs.  (1),  (3),  (16),  (18)  and  (22)  constitute  at  each  B  a  set  of 
five  algebraic  equations  in  the  unknowns  n,  u,  j, ,  n(o)  and  dVa 
(through  c* ) ,  provided  dVe  is  regarded  as  a  given  quantity  (here  taken 
to  be  independent  of  x),  and  provided  B(o)  (or  I)  and  V  are  both 
specified.  If  this  set  can  be  solved  for  each  value  of  B  between  B(o) 
and  a  value  to  be  assigned  to  the  channel  end  (we  will  take  this  to  be 
zero,  as  in  Ref.  (1)  ),  then  the  special  variable  x  can  be  calculated 
through  integration  of  Aapere’ ■  law: 


x 


dB 
Mo  3* 


(23) 


and,  in  particular,  the  full  channel  length  is 


More  rationally,  L  would  be  a  specified  channel  dimension,  and 
only  one  of  the  electrical  terminal  properties  (I  or  V)  would  be  given. 
This  will  require  a  global  iteration  in  which  terainal  values  are  tried 


until  the  proper  length  is  achieved. 

3.  Method  of  Solution 
Froa  continuity  and  aaaentua: 

U  =  f  1(0)*  -  B*  ]  (25) 

and 

n  =  -A-  (26) 

ai  Au 

These  expressions  can  be  substituted  into  Bq.  (18),  resulting  in  a 
relationship  linking  J,  with  B  and  4Vt  (through  c*)»  Algebraic 
rearrangeaent  yields  a  cubic  equation  for  j,  ,  which  can  be  coapacted 
to  the  fora 


(27) 

(28) 
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b  =  -  (1  +  2^+12^  + 

_  .  12  kT  _ 1 

~  b»  [B(0)«-  B*  ] 

and,  after  solution  for  {, 


(29) 


(30) 


J,  =  C-*B*  «  -  1  ♦  |  n)  (31) 

ht  a1 

Since  q  >  0,  the  cubic  discriainant  y  +  ^  can  be  shown  froa  (28) 
and  (29)  to  be  positive,  which  indicates  a  single  real  root  for  Bq. 
(27).  This  root  is 


+  [ 


b 

2 


✓ 


b1  a1 
4  +  27 


]i/i 


(32) 


The  algoritha  for  solution  now  involves  the  following  steps: 

(1)  Prescribe  I,  H,  L,  d,  a,  4VC 

(2)  Guess  V 

(3)  Select  values  of  B  between  B(o)  =  Mo I/H  and  o 

(4)  At  each  B,  guess  4Va ,  calculate  c* 

(5)  Coapute  q,  a,  b  and  (,  leading  to  j, 

(6)  Coapute  V  froa  Bq.  (21).  If  not  right,  aodify  AV  ,  return  to 
(4) 

(7)  When  converged,  update  B  by  B  •*  B  -  IB,  advance  x  by 


(8)  Continue  to  B  =  0,  obtain  L  =  x  (B  =  0) 

(9)  If  L  not  right,  aodify  V,  return  to  (2). 

(10)  Iterate  to  convergence. 
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4  Conditions  for  Onaet 

As  I  is  increased  at  a  fixed  m,  it  will  be  found  that  the 
increased  cross-force  j,  B  leads  to  a  progressive  reduction  of  density 
at  the  anode  (n(o)),  requiring  a  reduction  of  the  retarding  potential 
drop  (a  less  negative  AV, )  in  order  to  satisfy  the  increased  jv  demands 
(See  Bq.  (14)).  An  alternative  is  for  itself  to  be  reduced,  but, 
clearly,  this  cannot  happen  over  the  whole  channel  length,  since  the 
net  I  is  increased;  it  say  happen,  however,  over  the  central  portion  of 
the  channel,  where  u  B  nearly  cancels  the  whole  core  field,  but  u  B  is 
snail  at  the  channel  ends,  where  either  u  or  B  are  saall.  As  a 
consequence,  we  obtain  a  progressive  redistribution  of  current  towards 
the  two  ends  of  the  channel,  with  dVt  simultaneously  becoming  zero  or 
positive  over  the  center,  and  negative  only  over  smaller  and  smaller 
sections  near  the  ends.  These  sections  have  axial  thicknesses  of  the 
order  of  the  channel  length  divided  by  the  magnetic  Reynolds  number 
based  on  length,  and  are  therefore  rather  narrow.  Thus,  the  eventual 
consequence  is  an  inability  for  the  channel  as  a  whole  to  deliver  the 
prescribed  current,  a  steepening  of  V(I)  and  a  saturation  of  I. 
Clearly,  any  further  increase  in  current  must  imply  a  radically 
different  conduction  mechanism,  most  likely  in  the  form  of  highly 
erosive  vacuum  arcs  on  the  most  loaded  parts  of  the  anode,  so  as  to 
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provide  extra  aase  through  which  the  current  aay  channel  i'taelf .  This 
is  the  "onset”  condition. 

In  teras  of  our  aodel,  we  can  obtain  (following  Ref.  (1)  )  a 
simple  condition  for  onset  by  noting  that  n(o)  becoaes  then  ouch 
smaller  than  n  or  n(d).  In  Bq.  (17)  this  Beans  that  the  first  tera 
inside  the  root  is  negligible,  and  then  Bq.  (18)  siaplifies  to 


□  = 


oB1  j,  d 
ekT 


Solving  for  j,  and  using  (25)  and  (26), 


_  «  ekTtti  1 

-  9  OB'#*! 


B*  (B(o)*-  B* ) 


(33) 


(34) 


This  explicit  relationship  between  j,  and  B  can  then  be  inserted 
into  Bq.  (24)  and  integrated.  Rearranging  the  result  we  obtain: 


11  1  945  ekT  LH* 

i1  ~  8  om£  d* 


(35) 


For  fixed  channel  dimensions,  and  ignoring  the  relatively  weak 
variations  of  T  and  o,  this  relationship  has  the  fora 

Ii .  t i  - 

— t —  y  B|  :  const 


which  is  very  close  to  the  experimentally  observed  onset  condition 
I*  — 

j"  =  const,  (and,  incidentally,  also  to  the  similar  expression 
obtained  by  equating  the  exit  velocity  to  the  Alfven  critical 
ionization  velocity). 

As  an  example,  we  have  used  values  close  to  those  that  would  apply 
to  a  standard  Princeton  MPD  thruster  (with  some  uncertainty  about  the 
appropriate  channel  width  d  to  be  used  for  a  coaxial  geometry) : 

H  =  2n  x  5  cm  L  =  15  cm  d  =  H,-  ^  =  8  cs 

T  =  30,000  K  a  =  8000  O'1/*  (Argon). 

For  m  =  6  g/sec,  we  then  calculate  an  "onset"  current  of  14,900 
Amp.  This  is  somewhat  lower,  but  of  the  same  order  as  those  measured 
experimentally  (20-25  KA). 

In  performing  numerical  calculations  using  the  more  general  (below 
onset)  conditions,  we  have  indeed  observed  current  saturation  at  a  I 
value  corresponding  very  closely  with  that  given  by  Eq.  (35).  These 
calculations  are  reported  in  the  next  section. 

5-  Sample  Results  Below  Onset 

To  illustrate  the  features  of  MPD  flows,  we  have  obtained 
numerical  results  from  the  theory  explained  in  Sections  2  and  3  for  two 
examples,  defined  as  follows: 


Case  A:  Argon  gas, 


m  =  6  g/sec 
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T  =  20,000  K  a  =  8000  ST1/* 

H  =  31.4  ca  d  =  8  ca  L  =  11.2  ca 

Came  B:  Argon  Gas  m  =  6  g/mec 

T  =  20,000  K  a  =  5000  Q-'/n 

H  =  31.4  csa  d  =  5  ca  L  =  27.8  ca 

The  calculated  onset  current,  froa  Eq.  (35)  is  13.6  KA  in  Case  A 
and  23.03  KA  in  Case  B,  the  difference  being  aoetly  due  to  the  narrower 
gap  (d)  in  Case  B. 

For  Case  A,  Fig.  3  shows  the  axial  distribution  of  the  negative 
anode  drops  at  three  current  levels,  one  of  which  is  very  near  onset. 
Notice  that  for  the  given  conditions  the  anode  "floating  potential" 
(i.e. ,  the  AVa  for  which  j,  =  0  according  to  Eq.  (14)  is  -8.05  Volt. 
We  see  in  fig.  3  that  at  10  KA,  and  even  at  12.5  KA,  AVa  stays  within  1 

volt  of  this  level,  whereas,  as  we  coae  very  close  to  the  onset 

current,  AVa  rapidly  approaches  zero  or  even  reverses  sign  (although 
this  happens  so  near  onset  that  calculations  are  difficult  in  that 
case).  The  saae  trend  is  even  aore  clearly  apparent  for  Case  B,  in 
Fig.  6.  One  can  also  notice  the  survival  of  the  large  negative  anode 
drop  (even  near  onset)  in  the  channel-end  regions.  This  is  where  the 
back  eaf  uxB  is  saall,  and  the  anode  still  has  to  keep  up  its  barrier 
to  counteract  nearly  the  full  interelectrode  potential.  The  Magnetic 
Reynolds  nuabers  based  of  the  width  of  these  regions  (about  lea  in 
either  case)  is  about  0.5,  as  was  to  be  expected.  The  length-based 

aagnetic  Reynolds  nuabers  are  7  and  14  for  Cases  A  and  B  respectively. 
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center  region  of  the  channel  is  relatively  unaffected  aa  the  net 
current  increases,  since  the  corresponding  channel  voltage  increase  is 
being  compensated  by  a  similar  increase  in  AVt .  On  the  other  hand,  the 
end  regions  are  strongly  affected,  and,  to  a  first  approximation  can  be 
seen  to  collect  all  of  the  extra  current  as  one  approaches  onset.  The 
same  trends  can  be  seen  in  Fig.  7  for  Case  B.  Of  course,  the  current 
peak  at  the  downstream  end  is  to  be  interpreted  as  a  "bottled-up” 
reflection  of  what  normally  would  be  a  two-dimensional  current 
distribution  extending  into  the  plume;  we  have  forced  it  back  into  the 

channel  by  imposing  B  =  0  at  the  exit  plane. 

Figs  5  and  8  (Cases  A  and  B  respectively)  summarize  the  salient 
trends  versus  total  current.  Shown  are: 

(a)  The  total  transverse  voltage.  Notice  the  similar  shapes  in 
both  cases,  with  a  sharp  upturn  as  the  onset  current  I*  is  approached. 

(b)  The  negative  anode  drop  (at  aid-channel),  showing  the 
constant  value,  close  to  floating  potential,  up  to  about  90k  of  I"  , 
followed  by  the  rapid  drop  (or  increase  in  +  dV, ) .  Even  if  the 
absolute  value  of  the  computed  dVt  may  be  unrealistic,  due  to  the 
neglect  of  near-wall  excess  Ohmic  losses,  this  trend  should  be 
observable  at  currents  approaching  1* . 

(c)  The  plasma  densities  at  the  anode  and  cathode  walls  (at 
mid-channel).  Here  one  sees,  first  of  all,  the  general  reduction  of  n 
aa  current  increases,  which  is  due  to  the  higher  velocity  produced  and 
the  constancy  of  m  and  A.  But,  more  interestingly,  Figs.  5  and  8  show 


the  dramatic  departure  of  i^IODI  fro*  Dcathode  a®  I  approaches  I1,  and, 
in  fact,  the  approach  of  B0 D B  to  zero  at  I*.  Due  to  this  radical 
redistribution  of  *ass,  the  near-cathode  density  actually  stops  falling 
with  current  in  the  regime  near  I* . 

6.  Conclusions 

The  results  presented  above  indicate  several  interesting  facts: 

(a)  The  syaptoas  of  "onset"  are  largely  absent  up  to  some  90*  of 
the  onset  current,  which  shows  the  desirability  of  operation  fairly 
close  to  I*  . 

(b)  The  anode  drop  should  remain  nearly  constant  along  the 
channel  length,  except  at  currents  near  I"  .  This  indicates  that  si*ple 
models  based  on  an  "equipotential"  anode  should  be  acceptable  for 
I  S  0.9  P  . 

(c)  The  rapid  increase  of  4Vt  predicted  near  onset  should  be 
observable  by  probing  the  plasaa  potential  near  the  anode  surface. 
This  should  constitute  a  strong  verification  of  the  line  of  modeling 
taken  here  and  in  Ref.  1. 

(d)  A  similar  consent  applies  to  the  near^-anode  density,  whose 
rapid  decrease  with  current  near  onset  should  be  detectable 
spectroscopically  or  by  means  of  probes.  This  has,  in  fact,  been  done, 
by  Rudolph  (Ref.  2). 
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Fig.  6.  ARGON,  in  =  6  g/sec  T  *  20,000  K  a  =  5000  Q  -1/m 
H  =  31.4  cm  d  =  5  cm  (I  *23, 031  A) 
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In  this  Appendix,  a  two-dimensional  induction  equation  for  MPD 
channels  will  be  developed  and  used  in  a  local  description  of  the 
thruster  entrance. 

For  this  work,  the  plasma  is  modeled  as  a  steady,  fully  ionized, 
single  temperature  gas  characterized  by  a  constant  electrical 
conductivity,  a,  and  a  Hall  parameter,  0,  which  is  proportional  to  the 
magnetic  induction  B  and  inversely  proportional  to  the  gas  density,  p: 

a  a  Constant  ;  0  s:  c  jj  (c  =  SSL  )  (1) 

where  mi  is  the  ionic  mass  and  e  is  the  fundamental  charge.  The  plasma 
is  assumed  to  follow  the  ideal  gas  law,  P  =  2nj  kT,  where  2n(  =  n,  +  n, 
is  the  total  number  density,  k  is  the  Boltzmann  constant,  and  T  is  the 
plasma  temperature.  Experimental  data  *  and  theory  1  have  shown  that 
the  electron  temperature  varies  slowly  throughout  most  of  the  discharge 
except  in  the  initial  ionization  region.  Therefore  the  discharge  is 
assumed  to  be  isothermal.  This  approximation  is  also  justifiable  due 
to  the  large  electron  thermal  conductivity  and  the  excess  available 
ionization  and  excitation  energy  found  in  the  MPD,  but  it  is  ultimately 
just  a  convenient  device  to  facilitate  the  analysis. 

The  generalized  Ohm’s  law  can  be  written  as 

o  (S'  +  ;&-)  =  3  +  J  +  £  ;  8’  =  B  +  u  x  B 


(2) 
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where  Jf  =  fi  B/B  and  fi  is  the  Hall  parameter,  G  is  the  aean  aass 
velocity  vector,  B  is  the  aagnetic  induction  vector,  B  is  the  electric 
field  and  J  is  the  current  density. 

We  will  also  sake  use  of  the  conservation  equations  for  aass  and 
aoaentua,  as  well  as  of  Aapere’s  law  and  B  -  field  irrotationality. 
For  a  two-diaens ional  geoaetry,  with  OZ  perpendicular  to  the  plane  of 
the  flow  and  OX  along  the  channel  axial  direction,  these  are 

vpG  =  0  (3) 

p(G-v)G  =  -  v  (p  +  2^-)  (4) 

3  =  v  x  L  ;  B  =  B  I,  (5) 

V  x  B  =  0  (6) 

In  Eq.  (4) ,  the  Lorentz  force  3  x  B  has  been  converted  to  a 
aagnetic  pressure  gradient  -  v(B*/2ua).  An  entirely  siailar 

transforaation  gives  for  the  Hall  field  in  Bq.  (2). 

;  J  *  »  *  -  sr  »(sr>  <7> 


and  using  the  aoaentua  equation, 

h  3  x  3  =  dr  C  5  +  *(2n,kT)]  (7b) 

u  ezi| 

Substituting  this  into  Oha’s  law  (Bq.  (2)  ) 


o(8*  + 


)  =  3  +  zr  t  +  2v(n,  kT) ] 


(8) 


; mi  "■‘■I 
enj 


It  can  be  seen  that  the  effect  of  the  electron  pressure  gradient 
ten  in  Oha’s  law  is  to  cancel  1/2  of  the  siailar  ten  appearing  on  the 
right  hand  side  as  a  consequence  of  the  Hall  field  (sore  generally,  in 
a  partially  ionized  gas,  it  would  cancel  the  vP,  and  leave  the 
v(Pt  +  P, )  part).  We  can  rewrite  Eq.  (8)  as 


3  =  o[8‘  -  (a-y)  J  -  -  7  In  n,  ] 

6  6 


(9) 


where  the  (-)  sign  in  the  vln  nj  ten  is  to  be  carefully  noted. 

It  is  useful  to  think  of  the  whole  bracketed  ten  in  Eq.  (9)  as 
the  effective  electric  field.  This  field  contains  parts  which  are 
irrotational,  plus  others  which  are  of  the  fon  u  x  (  ).  We  can  obtain 
this  separation  by  noting  that  (u-v)u  =  v(u*/2)  +  5  x  u,  where  w  is  the 
flow  vorticity,  S  =  V  x  u.  Using  also  8'  =-v$+uxB  ($=  electric 
potential),  we  obtain 

3  =  o  [-  +  ~lnnt)  +  ilx  (8  +  ^-3))  (10) 

Notice  here  the  grouping  together  of  8  and  jj1-  w  ,  both  of  which 
are  vectors  along  OZ.  It  is  interesting  that  the  ion  gyro  frequency  Q 
is  eB/^  ,  so  that 
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H  +  a.5  =  SL(i}  +  5)  (11) 

6  " 

Bq.  (10)  is  a  convenient  point  of  departure  for  obtaining  the 
induction  equation,  for  this  purpose  we  take  the  curl  of  both  aides 
and  use 

n]  =  7»xvxi:-  -  v*B 

0o  Po 

and 

v  x  (u  x  B)  =  S(yfl)  +  (B*v)u  -  (u*v)B  -  B(vu) 

=  -(u*v)B  +  uVp  =  -  p(G‘V) (— ) 

P  p 

with  an  entirely  similar  expression  for  v  x  (u  x  5) .  Here  we  have  used 
the  two-dimensional  ity  to  cancel  the  term  (B-v)u.  With  these 
substituRLons,  the  induction  equation  is  found  as 


hr  =  p(«-v)( 

<Vg 


B 


+•  5 


(12) 


At  high  Magnetic  Reynolds  number  (R*  =  pgoul),  Bq.  (12)  indicates 
that  the  quantity  —  (B  +  - )  is  convective  (conserved  along  a 

r  © 

streamline).  Alternatively,  we  may  say  that  ~ w  is  convective. 

P 

Under  some  conditions  —  itself  may  also  be  convective  or  constant. 

P 

D 

x.e.,  when  the  flow  is  uniform  in  a  cross-section;  in  such  cases,  -  is 

P 

a  convective  quantity  too,  and  this  property  will  be  exploited  in 
theoretical  work  now  in  preparation. 


To  conclude  this  section  we  note  that  an  alternative  fora  of  the 
induction  equation  can  be  obtained  by  straightforward  application  of 
the  curl  operator  to  Bq.  (2),  with  the  result,  in  two  dimensions, 


v*B  +  I 


£B  gg  _  *B 
dy  dx  dx  dy 


]  =  Mo<*P(«'*)(“) 


(13) 


where  the  bracketed  tera  is  equivalent  to  the  tera  ongp( u*v)(®-2)  in 
Bq.  (12),  and  indicates  that  whenever  B  and  0  (or  B  and  n( )  are 
functionally  related,  such  that  the  corresponding  gradients  are 
aligned,  we  will  then  have  pure  convection  of  vorticity  (convective  ^). 

t 

2.  The  Inlet  Region 

The  propellant  is  injected  at  low  speed  through  the  backplate  of 
an  MPD  thruster.  As  long  as  (a)  the  aoaentua  tera  in  Bq.  (4)  is  saall 
coapared  to  either  the  Lorentz  force  or  the  pressure  gradient,  and  (b) 
the  back  eaf,  uB,  reaains  saall  compared  to  the  applied  field,  then  the 

general  Oha’s  law  (Bq.  (10))  reduces  to 

i 

i 

j  3  -  ~  °  *(♦  +  j—  In  nt )  (14) 


while  the  aoaentua  equation  reduces,  to  first  order,  to 


2n{  kT  + 


B* 

2Po 


=  const. 


(15) 


and  (as  can  be  seen  by  taking  the  curl  of  (14)),  the  induction  equation 


reduces  to 


v*B  =  0 


.16) 


Note  that  Eq.  (16)  does  not  necessarily  imply  a  constant  Hall 
parameter  0.  The  reason  the  bracketed  term  in  the  aore  general  fora 
Bq.  (13)  is  now  absent  is  that,  to  order  %  ,  B  and  n(  are  indeed 

functionally  related  (through  the  hydrostatic  balance  (15)  ).  It  is  an 
interesting  fact  that  the  very  siaple  law  (16)  for  B  is  found  to  apply 
even  in  the  presence  of  potentially  large  Hall  paraaeter  gradients. 

Using  Bq.  (7b),  Oha’s  law  (14)  becoaes 

3  +  3  x  (|)  2  oB  (17) 

indicating  an  effective  Hall  parameter  reduced  by  1/2  due  to  the 
electron  diffusion  effect. 

There  are  several  issues  we  would  like  to  explore  in  connection 
with  this  inlet  region: 

(1)  The  corner  singularites  in  B,  i.e. ,  current  concentrations, 
due  to  the  action  of  the  Hall  effect.  This  is  siailar  to  the  known 
current  concentration  effect  on  the  downstreaa  edge  of  cathode  segments 
in  MHD  channels,  except  that  in  the  case  of  an  accelerator  (uB  <  By ) 
the  concentration  will  be  at  the  upstream  cathode  edge. 

(2)  The  redistribution  of  mass  in  response  to  the  lateral  3  x  8 
force  as  it  develops  away  froa  the  back  plate  (at  the  back  plate 


itself,  3  x  8  ii  purely  axial,  though). 


(3)  The  ionization  and  heating  effects  associated  with  flow 


entrance. 


Bffects  (2)  and  (3)  are  for  now  deferred,  and  we  will  concentrate 


on  effect  (1)  above.  For  this  purpose,  Laplace's  equation  (Bq.  16) 


must  be  solved  with  appropriate  boundary  conditions.  In  formulating 


these  boundary  conditions,  we  will  assume  nonzero,  but  constant  Hall 


parameter  0.  This  is  permissible  because  only  the  cathode  corner  is  of 


interest,  and  the  fact  that  the  value  of  0  or  the  anode  wall  may  be 


different  (higher)  is  of  little  significance,  because  very  little 


current  will  in  fact  go  to  the  upstream  anode  corner.  In  what  follows 


we  will  implicitly  assume  a  constant  Hall  parameter  0  of  the  order  of 


(1/2)  the  value  at  the  cathode  corner. 


For  solution,  this  section  of  the  MPD  will  be  modeled  as  a 


constant  area  channel  at  a  constant  potential  (Fig.  1).  In  the  absence 


of  strongly  variable  sheaths,  the  tangential  electric  field  at  the 


anode  and  the  cathode  is  zero.  Additionally,  at  axial  distances  far 


fraai  the  backplate,  the  axial  electic  field  goes  to  zero  and  the 


transverse  electric  field  goes  to  V/h.  These  assumptions  lead  to  the 


boundary  conditions  for  the  channel 


at  (x,y)  =  (0,y) 


.  _  1  dB  _  n 

j*  •  1  ^  •  0 

1  mo  <*y 


at  (x,y)  =  (x,0),(x,h) 


dB  dB 
dy  "  p  dx 


at  (x,y)  =  (x  «*,y) 


where  uB  has  been  retained  despite  the  inconsistency  with  the  basic 
assumptions  explained  above,  since  this  solution  is  then  valid  for  a 
constant-velocity,  constant  fi  channel. 

This  is  a  well  posed  elliptic  formulation  that  can  be  readily 
solved.  In  the  complex  representation,  the  complex  variables  are 
expressed  as 

z  =  x+iy 

B  =  B  +  it  (19) 

The  complex  magnetic  field  and  the  magnetic  field  in  real  coordinates 
can  be  connected  through  the  definition  of  the  complex  derivative 

A  A 

dB  dB  .  dB  dB  dw 
d z  ~  dx  1  dy  ~  dw  dz 

Bquation(20)  shows  that  the  magnetic  field  gradient  consists  of  two 
parts,  one  relating  the  original  coordinate  system  (z)  to  a  new  mapped 
system  (w)  and  another  describing  the  variation  of  the  complex  field  in 
the  mapped  system. 

The  mapping  part  of  the  solution  can  be  facilitated  by  the 
Schwartz-Christoffel  Transformation.  This  provides  a  simple  mechanism 
for  converting  the  semi-infinite  region  defined  in  the  real  (z-plane) 
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to  the  upper  half  of  the  coaplex  (w-plane)  bounded  by  the  real  axis 


(see  Fig.  2) 
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Fig.  1.  Channel  Geometry 
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Fig.  2.  Coordinate  Transformation 
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This  transformation  is  accomplished  using  the  function 


w  =  cosh  (g5)  (21) 

Construction  of  an  analytic  function  relating  the  complex  magnetic 
field  to  the  coordinate  system  in  the  w-plane  requires  mapping  the 
boundary  conditions  from  the  z-plane  to  the  w-plane.  In  the  2-plane, 

Jn 

the  boundary  conditions  on  ^  can  be  obtained  readily  from  Bqa.  (18) 
and  (20): 

Anode  and  Cathode:  Arg  (~)  *  Arg(l  -  ifi)  -  v  -  tan' 1 0 

A 

Jn 

Backplate:  Arg  (gg)  *  v 


the  mapping  (21)  subtracts  o  on  the  anode,  on  the  backwall  and  ir 

A 

Jn 

on  the  cathode,  leaving  the  -jjj  arguments  shown  in  Fig.  3,  which  can  be 
mathematically  expressed  as 


along  A'B' 
along  D'B1 
along  B' O' 


arg(^-)  =  n  ~  tan' 1  (3 


arg(g-)  =  -  tan" 1 0  (22) 


If  the  analytic  function 


XT  =  K(w  -  l)“(w  ♦  1)T 


(23) 


is  selected  with  K  being  a  c  cap  lex  constant,  the  coefficients  a  and  t 

A 

and  the  argument  of  K  can  be  related  in  order  to  satisfy  (22),  and  are 
and  found  to  be  equal  to 


Using  Bq.  20  it  can  then  be  shown  that 


cut  n 


(27) 


Solving  for  the  real  and  imaginary  parts  of  Bq.  (27),  and  referring  to 
Bq.  (20),  the  partial  derivatives  of  the  aagnetic  induction  can  be 
derived  as  follows 


2®  =  _  a.Ki 

dy  h,E| 


,Jh c 


r cosh  ( )  —cos  ( jp" ) 


[cosh(jp)+cos(|j*)J 


sin[(tan'l0)(l 


sin(jp^)  l 


L.inh(p) 


)] 


(28) 

A 

For  x  -*  ®,  the  constant  |K|  can  be  calculated  froai  the  boundary 
conditions,  Bq.  (18),  and  is  found  to  be 


|K|  =  ^(1  +  /*)-*'*(£  -  uB) 


(29) 


The  above  results  show  that  a  closed  fore  solution  for  B  say  not  be 
possible;  however,  these  results  are  exact  expressions  for  the 
components  of  the  current  density  vector.  At  the  insulator  backplate. 


Ji  =  0 


^2tan~ 1 

j,  =  o(l  +  0*  )-"*(£  -  uB)[tan(g*)J  n  (30) 

( 

The  distribution  of  j,  along  the  backplate  can  be  seen  in  Fig.  4.  The 
anticipated  singular  region  is  seen  at  the  cathode  root.  Similarly, 
along  the  anode  (y=0) 


Ji 


=  -0o(l  ♦  P  )’* (^  ~  uB) [tanh(|^)J 


,2tan~ *0. 
1  U  } 


Jv 


V  ffv 

=  0(1  +  0*)-l(r  -  uB)  [ tanh(sr) ] 


2tan~ 1 0. 
(  It  ' 


(31) 
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and  along  the  cathode  (y=h) 

^2tan~ 1 0^ 

j,  =  -fio(  1  +  -  uB)[coth(§*)]  n 

^2tan-1d^ 

j,  =  o(l  +  *  )- »  (£  -  ufl)  [coth(§*)  ]  n  (32) 


The  diatribution  of  current  at  the  anode  and  the  cathode  is  diagraaaed 
in  Figs.  5  and  6. 

It  has  been  shown  that  a  very  sisple  aodel  of  the  entrance  region 
of  the  MPD  can  be  constructed  that  has  a  lisited  range  of  validity  but 
deannstrates  the  singular  nature  of  the  channel  current  distribution. 


i 

i 

i 

i 
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Cothode 


Fig.  4. 

Distribution  of  transverse 
current  density  on  the  face 
of  the  Backnlate 
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Fig.  5. 

Distribution  of  transverse 
current  density  along  the 
anode  wall 


Cathode 


Fig.  6. 

Same  as  Fig.  7,  along  the 
cathode  wall 
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Appendix  3 


Simultaneous  Hall  and  Magnetic  Convection  Effects 

Consider  a  rectangular  channel  (Fig.  1)  of  length  L,  height  W 
and  depth  H,  with  infinitely  segmented  electrode  walls  and  bounded  by 
electrical  insulators  at 


PiO.  i  .  Aet-r*  AA 

CuArfrtf*-  SC  *AT  I  t. 


its  ends.  A  fluid  of  constant  conductivity  o  and  Hall  parameter  3 
is  assumed  to  flow  in  the  channel  at  a  constant  velocity  u.  We  define 
as  a  reference  non-dimensional  quantity  the  Magnetic  Reynolds  number 


R  »  yQa  u  L 


(1) 


where  ■  Ait  x  10  ^  (MKS  units)  We  also  normalize  distance  as 


w 

w  “  — 
L 


(2) 


If  we  impose  the  design  condition  that  the  current  vector 
j  be  everywhere  transverse  (j  =  0),  it  can  be  shown  that  the 


J  component  and  the  potential  $  with  respect  to  the  cathode  tip  are 
given  by 


K  *  a  Bo  u  T“ , 
3  e  -1 


♦  RtT-n+|a-«-R(1  -O) 

1  -e 


where  Bq  is  the  magnetic  field  at  x  ■  0,  related  to  the  total 


current  I  by 


I-S  B 

wo  ° 


Eq.  3  shows  that,  if  R  >>  1,  most  current  will  concentrate  in 
the  last  1/R  fraction  of  the  channel  length. 


/S-*i 


f-e" 


4  % 


voLTA&ex  r»R 
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If  there  ie  substantial  overlap  between  the  anode  and  cathode 
voltages  (Fig.  2),  i.e..  If 


6  » 


0 

1  -e 


u  W 

'-IT 


(6) 


then  the  required  voltage  distribution  could  be  approximated  using 
a  single  external  supply  by  externally  interconnecting  points  on 
the  anode  and  cathode  which  must  be  at  equal  potentials.  This  is 


schematized  in  Figs.  3  and  4, 


FiG.  3  u- 

2a~1c  ,  e-rc.  be 


fit  4  (.ftod-tONfjecn  o*/ 

ScH**jiri£  rox  Fi6.  t\3 


where  a  six-segment  construction 
would  result.  The  dotted 
lines  in  Fig.  4  Indicate 
external  shorting  wires. 

Notice  the  total  channel 
current  is  I  •  6  I  • 

CAL  • 

where  I  _  is  the  power 
ext 

supply  current,  while  the 

lc, 

external  voltage  V£x  is 

6(V  +  AV  +  AV  ),  where  V 

y  a  c  y 

is  the  channel  voltage,  and 

AV  and  AV  are  the  near- 
a  c 

electrode  voltage  drops.  The 
fractional  voltage  drop  is 
the  same  as  in  a  direct  connec¬ 
tion,  since  both  the  plasma  and 
and  the  electrode  drops  are 
multiplied  six  times. 
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Scaling  and  Numerical  Estimates. 

The  Hall  parameter  $  can  be  estimated  as  follows: 


8 


U  a, 

o  i  _  ..I 

- a  u  W  — 

e  • 

m 


(7) 


where  full  ionization  has  been  assumed,  and  Bq  has  been  used 
for  B,  since  the  magnetic  field  tends  to  fall  mostly  near  the 
exit.  Condition  (  6)  becomes 


I 

e 

m 


>> 


(e/n^) 


1  -e 


-R 


(8) 


— R 

For  Argon,  assuming  e  <<  1, 


this  is 


>>2  4  — — - 

(g/sec) 


(9) 


This  can  be  compared  to  an  "onset"  criterion;  for  Argon,  an 
approximate  condition  for  stability  is 


-  <  82  7^ 

•  'V'  (g/sec) 

m 


(10) 


Assuming  0  •  4000  mho/m,  u  ■  6000  m/sec,  W  *  2  cm,  Eq.(  7) 


gives 


6  *  0.25  7 

m 


(KA) 

(g/sec) 


(ID 
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Fig.  5  shows  conditions  (  9)  and  (  10)  together  with  some 

constsnt-6  lines.  It  can  be  seen  from  this  figure  that  sufficient 

overlap  exists  for  m  below  some  2  g/sec,  even  if  the  ordinary  onset 

condition  is  not  exceeded  (of  course,  segmentation  is  intended  as  a 

means  to  exceed  that  condition).  Incidentally,  for  the  same  conditions, 

— R 

R  -  30  L(m.),  so  that  for  L  >  7  cm  the  assumption  e  «  1  is  justified. 

This  analysis  shows  that  axial  current  can  be  suppressed  by  appro¬ 
priate  voltage  distributions  on  the  electrode  walls,  and  that  such 
distributions  can  be  approximated  by  cross-connected  segmentation, 
using  still  a  single  power  supply.  It  also  shows,  however,  that  the 
bunching  up  of  current  at  the  channel  exit  will  persist  if  R  >>  1; 
for  an  open  channel  end,  this  would  manifest  itself  as  current  lines 
"ballooning"  into  the  plume. 
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CHAPTER  I  INTRODUCTION 


High  specific  impulse  with  moderate  thrust  density  and  relative  operational 
simplicity  are  the  main  features  of  the  quasi -steady  electromagnetic  plasma  accelerators 
(Magnetoplasmadynamics  MPD  Thrusters).  In  previous  research,  various  schemes  of 
application  of  external  magnetic  field  and  geometric  variations  were  tested.  [1]  One 
particular  concept  which  uses  coaxial  geometry,  high  current,  flow  rate  in  millisecond 
pulses  and  no  external  magnetic  field  [2.3],  however,  is  the  prime  interest  in  this 
country.  Considerable  progress  has  been  made  in  the  design,  operation  and 
of  these  thrusters.  Still,  two  areas  must  be  progressed  and  developed  in  order  to 
bring  these  features  of  the  thruster  into  functional  reality.  They  are  (i)  to  increase 
efficiency  [4]  and  (ii)  to  increase  lifetime  [S]  of  the  accelerator.  In  the  last 
many  test  data  regarding  the  operation  limitations  of  the  accelerator  have  been 
accumulated.  Most  of  them  aie  on  die  self -field,  coaxial  type  [6,7]  mentioned  above. 
(Figure  1.1,  1.2  &  1.3)  Important  efforts  have  also  been  made  in  the  research  of  its 
lifetime.  Not  much  data,  however,  has  been  collected.  More  efforts  are  now 
underway  [8] . 

Regarding  the  performance  limitation  analysis,  one-dimensional  flow  models  [9] 
have  provided  hints  about  improvements  of  geometrical  contouring,  but  the  essentially 
two-dimensional  phenomena  at  onset  of  fluctuations,  which  affects  the  attainable 
efficiency  [10] ,  cannot  be  addressed  by  these  analysis.  Some  theories  [11]  attempt  to 
correlate  the  onset  phenomena  by  invoking  the  "Alfven"  limiting  velocity  in  a 
dimensional  manner.  Highly  simplified  two-dimensional  theories  have  also  been 
developed,  [12]  but  geometrical  shape  of  the  thruster  still  remains  to  be  the  deciding 
factor. 
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A  physically  illuminating  analysis  by  Baksht,  et  al  [13]  postulated  that  the 
side  force,  resulting  from  axial  Hall  current,  reduces  the  near-anode  plasma  density. 
Then,  current  saturation  will  be  reached  when  the  whole  anode  has  positive  electric 
potential  relative  to  the  plasma.  Therefore,  other  phenomena  must  occur  if  higher 
current  is  demanded. 

Other  phenomena  accompanying  current  satuation  are  full  ionization  &  high 
Hall  parameter  near  the  anode,  and  partial  ionization  but  high  electron  &  neutral 
number  density  near  the  cathode.  These  microscopic  properties  suggest  that  the 
analysis  of  local  instability  (microscopic  instability)  might  provide  answers  for  the 
global  onset  situation.  Since  the  Hall  parameter  and  electron  density  play  a  major 
role  in  local  stability.  More  understanding  regarding  those  parameters  is  needed. 

Therefore,  this  thesis  will  discuss  the  local  instability  of  the  plasma  inside 
the  coaxial  accelerator.  In  particular,  a  small  disturbance's  propagation  through  the 
plasma  of  uniform  background  will  be  the  focus  of  our  attention.  Instabilities 
resulting  from  ionization  and  field  variations  will  also  be  discussed  in  the  later 
chapters. 


CHAPTER  II  PROPAGATION  OF  DISTURBANCES  THROUGH  THE  PLASMA 
WITH  PARTIAL  IONIZATION.  CONSTANT  TEMPERATURE  &  CONSTANT 
INDUCED  MAGNETIC  FIELD 

LI  INTRODUCTION 

Global  parameters  such  as  electron  temperature,  density  are  important  as 
premises  for  the  formulation  of  our  local  analysis  since  they  determine  the 
background  properties.  Before  starting,  we  must  have  some  knowledge  regarding  the 
global  situation. 

Baksht,  et  al  [13]  analyzed  the  onset  phenomena  for  a  simple  slender 
channel  geometry  by  assuming  constant  temperature  and  full  ionization  throughout  the 
plasma.  The  constant  temperature  approximation  is  substantiated  by  the  measurements 
made  by  K.  Rudolph  [10]  and  P.J.  Turchi  [14],  but  full  ionization  can  only  be 
reached  in  some  particular  regions  of  the  thruster. 

Based  on  other  measurements  regarding  the  electron  density  and  temperature, 
the  plasma  is  at  a  density  level  one  to  two  orders  of  magnitude  below  the  Saha 
equilibrium  density.*  Therefore,  bulk  recombination  is  very  unlikely  to  compete  with 
the  diffusion  of  charged  particles  to  the  wall  of  the  thruster,  (recombination  can 
only  occur  at  the  wall) 


See  appendix  I 


Radiation  recombination  is  certainly  present,  but  the  process  is  slow  in 
comparison  with  diffusion.  It  might  have  some  effects  on  the  local  stability,  but  in 
this  chapter,  our  prime  interest  is  to  set  up  a  simple  model,  so  we  can  investigate 
the  basic  response  to  small  disturbances  created  inside  the  plasma.  Also,  because  of 
its  comparison  with  particle  diffusion  rate,  we  can  safely  ignore  radiation  recombination 
here. 

Induced  magnetic  field  is  actually  varying  with  time  and  space,  as  stated 
by  Maxwell’s  equations.  These  variations  are  first  assumed  to  be  slow  and  negligible 
in  the  dimensions  of  interest,  so  it  is  assumed  here  that  magnetic  field  is  constant 
The  effect  of  varying  magnetic  field  will  be  discuss  in  the  next  chapter. 

So.  our  first  model  for  analysis  should  follow  the  guidelines  of  the 
classical  glow  discharge  theory  of  Schottky.  [15]  where  bulk  electron  collision 
ionization  is  mainly  balanced  by  ambipolar  diffusion.  The  net  gain  of  electron  in 
ionization  process  is  caused  only  by  disturbances.  In  the  global  energy  balance,  all 
energy  obtained  by  the  plasma  from  external  sources  (by  means  of  Joule  heating) 
should  be  balanced  by  losses  due  to  particle  collisions. 


2.2  CALCULATIONS 


i 


2.2.1  ASSUMPTIONS 

From  the  global  picture  of  the  plasma  mentioned  above,  we  can  extract 
some  assumptions,  and  proceed  with  our  analysis  accordingly. 

1.  Constant  temperature  throughout  the  plasma. 

2.  Plasma  is  partially  ionized  but  no  net  gain  of  electrons  by  the 
ionization  process.  (No  steady-state  ion  densities  time  variations) 

3.  Constant  magnetic  field,  (no  space  St  time  variations) 

4.  Positive  ion  and  neutral  possess  the  same  value  of  temperature, 
which  remains  constant  but  differs  from  the  electron  temperature. 

2.2.2  GOVERNING  EQUATIONS 

For  the  plasma  of  interest,  the  electrons  have  very  rapid  motion  relative  to 
the  ions  and  neutrals! heavy  species),  simply  because  the  mass  of  the  electron  is 
substantially  smaller  than  that  of  the  heavy  species.  Also,  the  energy  transfer  by 
electrons,  by  whatever  means,  is  faster  than  any  other  species.  So  it  is  appropriate  to 
consider  the  energy  balance  of  the  electrons  only. 

The  energy  input  to  the  electron  by  Joule  heating  (  J 1  /  <r~ )  is  balanced 
on  average  by  the  collisional  losses  between  the  electrons  and  ions,  (at  moderate 
ionization  fractions,  electron -ion  collisions  are  much  more  frequent  than  electron-neutral 
collisions)  However,  in  the  area  where  disturbances  occur,  the  excess  energy  given  by 


the  imbalance  process  will  fuel  the  ionization  process.  Thus,  the  energy  equation 
becomes. 


-  4r  -  3^^KU.-Tf)  (21> 

where  «  Q./ mc^e.  is  the  collision  frequency  between  the  electrons  and 

ions.  And  /At.  "  /?>/&  is  the  mobility.  3  **  the  Hall  parameter. 

Electron  temperature  is  much  higher  than  the  positive  ion  and  neutral 
temperature  because  of  its  rapid  random  motion,  so  Te-Tg  -Te 


After  rearranging  the  energy  equation,  it  yields. 


3 


«.V. 


(2-2) 


Also,  the  momentum  changes  of  the  electrons  are  governed  by  Ohm’s  Law. 


(electron  momentum  equation) 


the  electromagnetic  relations  of  the  charged  particles,  reduce  to  the  form. 


V  *  t  >  o  (2.41 

^  -  J  *  “  (2.5) 

2.13  CALCULATION  PROCESS 

Before  we  rearrange  the  equations,  let  us  set  up  the  coordinate  system. 
The  disturbances  are  assumed  to  be  of  a  plane  wave  form  propagating  at  any  random 


direction.  We  fix  our  coordinates  as  follow;  (  Figure  (2.1)  )  The  angle 


is  the 


angle  between  the  direction  of  the  propagation  of  the  disturbances,  and  the  steady 
state  (background)  current  flowing  direction. 

By  applying  the  two  Maxwell’s  equations  (2.4.  2.5).  we  deduced  that. 

Jx  ■  Constant 


Ey  -  Constant 


(there  is  no  variations  along  the  wave  front.  i.e  =  o  ) 


Let  us  denote  the  quantities  tn  the  background  state  by  the  bracket  <  >. 
For  example,  the  current  density  to  the  zero  th  order  will  be  represented  by  <  J  >. 


”  •  T  ^ 


T1  r  r‘  “  r 


i.»  *  *  r  i 


J-  M  * /  V 


From  the  assumptions  mentioned  above,  we  also  found  that  there  is  no 
background  time  variations.  i.e.  (  <+-€>')  ~  Since  electron  temperature  and 

magnetic  field  are  assumed  constant,  every  quantity  related  to  these  two  parameters 
carries  the  same  value  as  the  background  situation.  One  of  those  quantities  is 
conductivity  of  the  plasma.  In  the  regime  which  we  are  interested  in.  the  collision 
process  is  dominated  by  Cloumb  collisions.  The  conductivity  under  this  condition  is 
almost  only  a  function  of  electron  temperature,  as  appeared  in  the  Spitzer  formula. 


0s*  -  o  oj  5  2 


d0  =  4s!  j  Is_ 

J  "6 


-i  *  ’  i  A 

^ ) 


We  approximate  A.  as  a  constant  for  our  analysis;  taking  Ln  A.  «  6  , 


obtain 


q~  -  2.55  x  1<T3  Tcl5 
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The  disturbances  are  assumed  to  be  small  in  magnitude  in  comparison  with 
the  background  quantities.  We  therefore  can  write. 

J  ■  <  J  >  *  6  j 
(  J  has  two  components,  Jx  dk  Jy  ) 

{$  -  <  £  >  +  &(* 
fJk.  -  <  ♦  ^A'e- 

is  the  first  order  variation  (perturbed  variation),  and  so 

-i£-  «  x 

<e> 

After  all  these  substitutions  have  been  made  in  the  governing  equations 
(12.  13.  14  and  13).  and  the  zero'th  order  terms  cancelled  out.  we  obtain  the  first 
order  perturbation  equation 


&.(*)  - 


1  n 

-  3  ~ ^  L  ^  +  </2> - 


where  '3  -  k  Te  /  is  the  square  of  the  isothermal  speed  of  sound  in 

the  ionized  gas  (also  called  the  Bohn  velocity). 

Vj  is  the  voltage  (energy)  required  for  the  ionization  process  for  a 
particular  substance.  (Argon,  for  example) 
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2.3  RESULTS 


To  solve  the  differential  equation,  we  assume, 

— 

$£>  -  Z 


(2.8) 


where 


3  W  B 

Vc  (T- 


If  zA  is  positive,  the  disturbances  will  die  out  as  time  progremes,  in 
other  words,  it  will  not  propagate  through  the  whole  plasma,  and  so  the  system  will 
be  stable.  The  deciding  factor  in  ^  is 


=2  ♦  <fi  > 

The  variations  of  this  factor  due  to  the  changes  of  the  current  angle  and 
the  value  of  the  Hall  parameter  are  plotted  in  figure  (2.2)  and  figure  (2.3). 


2.4  CONCLUSIONS 


If  the  direction  of  the  wave  vector  (direction  of  the  propagation  of  the 

disturbances)  falls  between  the  direction  of  the  current  flow  and  electric  fieid. 

stability  occurs.  It  is  also  known  that  high  Hall  parameter  means  that  the  current 
flow  direction  is  further  deviated  from  the  electric  field  vector,  and  this  occurs  at 
relatively  low  electron  density  in  the  case  of  constant  induced  magnetic  field. 

In  figure  (2.4),  the  situation  falls  into  the  stability  category.  In  the  high 
Hall  parameter  region,  current  density  is  higher  than  the  two  adjacent  regions.  As  a 
result,  energy  deposition  by  means  of  Joule  heating  (  J2  /  <7-  )  will  increase,  and 
consequently  the  ionization  rate  will  go  up.  More  electrons  will  then  fill  the 

originally  low  density  regions,  thus,  the  discrepancy  will  ease.  Uniform  background 
can  then  be  preserved.  In  out  calculations,  the  angles  of  current  flow  which  fell 

into  this  regime  (0  to  90  :  180  to  270  )  did  exhibit  stability  at  any  value  of  the 
Hall  parameter. 

Conversely,  in  the  case  of  figure  (2.5),  Joule  beating  in  the  originally  low 
density  regions  will  decrease.  The  discrepancy  of  electron  density  will  then  be 
magnified,  and  so  the  uniform  background  would  be  destroyed.  As  proved  by  our 
analysis,  the  directions  of  current  flow  which  had  the  angle  between  90  to  180*. 
270  to  340  #  relative  to  the  wave  vector  were  the  unstable  region. 
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Another  conclusion  from  the  stability  condition  (equation  2.9) 

{  2  «•  <  p,  >  SIN(  2  )  }  >  0 

is  that  the  threshold  value  of  the  Hall  parameter  is  2  (corresponding  to  the  worst 
orientation.  £  -135  ).  It  serves  as  a  reference  point.  The  value  of  <£>  is 
considered  to  be  high  if  it  has  a  value  greater  than  two.  to  be  low  otherwise. 

In  the  "Barrel"  luminous  region  of  the  MPD  Arcjet  observed  by  K. 
Rudolph  [10]  at  the  onset  of  fluctuations,  the  Hall  parameter  did  exceed  2.  (it 
reached  10  in  that  region)  This  suggested  that  local  instability  could  have  an  effect 
on  the  phenomena  of  the  onset  of  fluctuations.  However,  more  refined  analysis  (to 
be  presented  later)  indicates  that  heat  conduction  should  provide  some  damping  in 
the  regime  of  interest. 


CHAPTER  III  PROPAGATION  OF  DISTURBANCES  WITH  PARTIAL 
IONIZATION,  CONSTANT  TEMPERATURE  BUT  VARYING  MAGNETIC  FIELD 


3.1  INTRODUCTION 

In  chapter  two,  we  proved  that  local  instabilities  could  contribute  some 
effects  to  the  global  onset  phenomena.  The  imbalance  between  the  Joule  heating  and 
the  collisional  losses  by  the  electrons  will  cause  instability.  But  there  are  other 
mechanisms  which  can  cause  instabilities  or  stabilize  the  existing  problem.  So  in  this 
chapter,  we  will  explore  the  effects  of  the  wavelength  dependence  in  wave 
propagation. 

We  neglected  the  variations  of  induced  magnetic  field  in  last  chapter  for 
simplicity.  But  in  some  cases,  their  associations  with  the  growth  of  instabilities  might 
be  important  Specifically,  a  space  varying  magnetic  field  might  not  be  able  to 
penetrate  through  the  whole  plasma,  since  the  "skin  effect"  will  act  as  a  damping 
force  against  the  penetration.  This  will  cause  a  wavelength  dependent  damping,  where 
certain  disturbances,  whose  wavelengths  are  of  the  order  of  some  characteristic  values 
associated  with  the  skin  depth  with  growth  lime,  will  die  out  preferentially.  Since 
the  normal  size  of  the  MPD  device’s  dimension  is  of  the  order  of  10~1  to  10~2 
meters,  the  disturbances  of  interest  to  us  must  be  of  this  order. 

As  stated  earlier,  the  magnetic  field  is  actually  varying  with  time  and 
space,  as  expressed  in  Maxwell’s  equations.  So  we  will  keep  every  assumption  the 
same  as  in  chapter  II  except  the  magnetic  field.  Thus  we  can  analyze  the  effects  of 
the  wavelength  damping. 


3.2  CALCULATIONS 


3.11  ASSUMPTIONS 

We  maintain  every  assumption  in  chapter  two  except  that  the  magnetic 
field  has  time  and  space  dependence. 

1.  Constant  temperature  throughout  the  plasma. 

1  Plasma  is  partially  ionized  but  no  net  gain  of  electrons  or 

positive  ions  by  the  ionization  process.  (No  time  variations  in 
ions  densities) 

3.  Ions  and  neutrals  possess  the  same  value  of  temperature,  which 
remains  constant. 

4.  The  magnetic  field  has  space  and  time  dependence. 

3.12  GOVERNING  EQUATIONS 

We  start  off  with  the  same  energy  equation  as  in  chapter  two.  where 
excess  energy  obtained  by  the  imbalance  between  Joule  heating  and  collisional  losses 
fuels  the  ionization  process. 


eVi 


dJU 


The  temperature  difference  between  the  electrons  and  heavy  species  is  still 


assumed  to  be  large  compared  to  the  gas  temperature.  Thus, 


Te  -  Tg  “  Te 


Substituting  the  proper  expressions  relating  to  n*  . 


Ue.  *nd  /< 


(Bohm  velocity)  is  the  square  of  the  isothermal  speed 

ionized  gases. 

The  electron  momentum  equation  is 


0“ 


c 


J 


J 


x 


/3 


With  proper  substitutions,  it  yields 


7- 


J 


j  <  (/'(.B) 


I 

I 

1 

» 


Maxwell’s  equations  which  govern  the  electromagnetic  effects 
charged  particles  are. 


- 


i 


.  we  get 


(3.2) 


of  sound  in 


(3.3) 


(3.4) 


between  the 


(3.5) 


^  •  B 


o 


(3.6) 


1  *  B 


(3.7) 


3.13  CALCULATION  PROCESS 


Since  our  interest  is  the  wavelength  dependence  of  the  propagation  of 
disturbances,  a  plane  wave  form  dependence  of  these  quantities  is  assumed. 


t  -  X  ■  1 

ou.  -  £_ 


(3.8) 


where  U.  is  constant  with  respect  to  time  and  space  coordinates, 
represents  the  quantity  of  interest. 

The  coordinate  system  is  the  same  as  in  chapter  II.  where  the  wave  ts 
propagating  along  the  x-axis.  The  angle  &  is  defined  as  the  angle  between  the 
background  current  flow  and  the  wave  vector(direction  of  propagation).  The  magnetic 
field  vector  will  be  perpendicular  to  the  x-y  plane.  (Figure  (2.1)  ) 
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Then, 


&  U. 


1  tO  u_ 


(3.9) 


dK 


fcU 


“*  ‘  &  uC 


(3.10) 


Small  disturbances  are  still  assumed  to  exist  Every  perturbed  quantity  a 
tiny  when  compared  to  the  background  state.  For  example: 


£> 


</'*> 


<.«*  I 


i.vi-  = 


4  la 
U 


As  in  last  chapter,  there  is  no  variations  along  the  wave  front  therefore. 
°  for  whatever  quantity  u  represents. 


Rearranging  the  Maxwell's  equation  (3.6.  3.7).  we  get 


J*  -  Constant 

Again.  every  quantity  in  the  background  state  is  constant 
(  3L  <ruv>  o  .y  -  o  .  etc)  except  the  magnetic  field,  as  appeared 


in  the  Maxwell’s  equation  (3.7) 


y  *  6  *  jj-  > 

The  magnetic  field's  background  state  variation  in  space  equals  to  the 
average  current  flow  density  at  that  state:  however,  this  needs  not  concern  us. 
because  gradients  of  <  B  >  do  not  appear  in  the  formulation. 

j*9  ^  3  > 

In  this  model  our  first  order  variables  are  £y Ut.  .  s£>  -  St*  -  Other 
variables  such  as  bjy  .  etc.,  are  just  dependent  variables  of  the  three  mentioned 
above. 

Before  we  proceed  any  further,  let  us  define  some  non-dimensional 
quantities  in  order  to  simplify  the  analysis. 


The  non-dimensional  frequency  will  be  defined. 


M 

where  £  can  be  interpreted  as  the  time  for  Ohmic  dissipation  to  generate  an 
energy  density  equal  to  the  magnetic  energy  density. 


The  non-dimensional  wave  number, 


-•n  <s? 
bx  47  <j  > 


K  <>  £> 
</> 


(3.12) 


(  /< 


where  7\  is  the  characteristic  macroscopic  variation  length  for  B,  such  that  a 
cyclic  change  of  B  with  as  the  wavelength  and  with  amplitude  <  B  >  would 

generate  currents  with  amplitude  <  J  >.  The  non-dimensional  wave  number  can 
also  be  interpreted  as  the  inverse  of  the  magnetic  Renold’s  number. 


V  v  \W 


Finally,  the  parameter  measuring  the  energy  invested  m  the  ionization 


process. 


<  e>  > 


(3.13) 


can  be  interpreted  as  the  ratio  of  ionization  potential  to  Ohmic  voltage. 


y|/<»is  the  magnetic  permeability  of  free  space.  A1**  4  1  1°“7 


SI  units. 


get. 


From  the  y -direction  component  of  the  electron  momentum  equation,  wt 
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>  > 
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V  >  <■/  "  " 


<3.141 
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find 


After  the  zero' th  order  terms  of  the  energy  equation  arc  cancelled  out.  *« 


-  *  l  *  <  /  /  t  ^  V 

/  4  A  , 


J  «\ 


(3  13) 
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% u  non-dimensional  quantity  for  . 

it  with  aquations  (3.8)  (3.9)  (3.10)  (3.11)  (3.13)  and  (3.14).  we  get. 


A  - 


v  - 


■i  <+•  >  >  'A 


♦  if) 


(3.16) 


f'c 


Subaututinc  ttoe  proper  expreMons  into  aquation  (3.7).  it  yields 


V 


t 


'  *  ?>‘  $v 


(3.17) 


v  •;  7  c~,>  -  «■  f  ?  =  o 


From  Maxwell's  aquation  (3.3). 


: 


<  .  -  5,  It  -  ^ 


(3.18) 
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Now,  the  final  form  of  our  governing  equations  are 


i 

\ 


where 


i 
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A  c  c .  j  ') 
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A  (  1.  1  )  - 


A  (  1.  2  )  - 


A  (  1,  3  )  - 


A  (  2.  1  )  - 


A  (  2.  2  )  - 


A  <  2.  3  )  - 


A  <  3.  1  >  - 
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A  (  3.  2  ) 


1  J  1  f  >' 


A  (  3.  3  )  -  —  c  -TL 

For  a  non  trivial  problem,  the  determinant  of  the  coefficient  must  equal 
to  zero.  In  the  determinant,  _n  is  treated  as  variable,  while  the  rest  is  considered 
as  input  parameter. 

Since  plane  wave  dependence  of  the  disturbances  is  assumed. 

feu.  =  *!'**-**> 

.  e  ‘  '^7  -  J  "lT  5 

(  'C.  .  L  are  constants.) 

We  notice  that  the  disturbances  will  exhibit  unstable  behavior  if  the 


imaginary  pvt  of  i".  is  neiative.  (  '4*  “  assumed  to  be  real  ) 


Expanding  the  determinant  into  a  quadratic  form,  we  get. 

sx  -  t  n(  +  t  -*•  l|»’; 


-  c  '+^2-[  j<p >e*js  -  * 


(3.19) 


Our  prime  interest  here  will  be  the  imaginary  part  of  the  unknown  XI  . 

Looking  into  the  equation  in  more  detail,  let  us  first  establish  the 
relationship  between  the  oscillating  part  (  Cl  real)  and  the  damping  part  (Cl 
imaginary).  By  separating  the  equation  into  real  and  imaginary  part  and  setting  them 
both  equal  to  zero,  we  find 
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i 

* 


f  XL* 


(3.20) 

f>  )  -  t* 


(3.21) 
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Rearranging  aquation  (3.21).  it  yields 


i</y> 
c  ~ 

(3.22) 

^  </3  >  '  '  i  ^ - 

j  ^  -  4>* 

We  focui  our  in  ter  eat  oo  the  worst  onenuuons  retarding  subility  (the 
magnitude  of  pL  tmlg  is  the  lowest),  so.  by  differentiating  the  equation  with  respect 
to  £  and  setting  the  derivative  .~lz)  -o.  we  find  the  angles  of  the 

worst  orientations  to  be  given  by 


+£> 


(3.23) 


We  further  specify  our  point  of  interest  to  be  the  threshold  situation 
where  the  plasma  is  in  neutral  subility  (zero  growth  rate  for  the  disturbances),  i.e 
,  -  0  So.  the  worst  orientations  and  the  relation  between  the  oscillating 
part  and  damping  part  are 


f  C*  r>  &  & 


(3.24) 
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—  4*  <8> 
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From  equation  (3.24).  we  see  that  {  COT  (  2  &■ )  }  tends  to  zero  both  at 
short  wavelengths  (  ^  >>  1)  and  long  wavelength  (  v|><<  i).  At  the  value  ^  - 
Jz  >  /  5',  (of  order  of  unity  in  cases  of  interest).  {  COT  (  2  9  )  l  peaks 
at 

[  COT  (  2  f  )  )m$x  -  r  j  /  8  -  0.2 

Since,  as  will  be  seen,  the  range  of  ^  in  the  cases  of  interest  is  from  about  7  to 
700.  it  can  be  concluded  that  {  COT  (  2  0  )  -  0  }  is  a  good  approximation  for 
maximum  growth.  This  gives  b  •  135  :  the  alternative  choice  $  -  45'  leads  to 
minimum  growth.  Details  of  the  variation  of  9  worjt  for  the  range  of  wavelengths 
from  1  mm  to  0.1  m  are  given  in  figure  (3.1). 

Now.  combining  equations,  we  get  a  relation  between  the  non-dimensional 
wave  number  and  the  Hall  parameter. 


We  can  see  that  at  both  limits  where  ^  is  zero  or  approaching  infinity 
(both  extremely  short  and  long  waves),  the  value  of  the  critical  Hall  parameter  is  2. 
The  long  wavelength  region  (  ^<<  0  )  is  one  where  the  disturbances  have 
wavelengths  larger  than  the  distance  required  for  the  current  to  build  up  the 

existing  background  magnetic  field.  Thus,  for  a  self-field  device,  this  region  is  of 
dubious  physical  significance.  The  opposite  limit  (  4*>>  1.  short  wavelength)  can  be 
interpreted  as  follow:  the  penetration  distance  in  a  time  ^  for  a  current 
disturbance  in  a  plasma  of  conductuvuty  >  is 

-  .  rr-’ 

i  ^  , - 

,  A'*  <f 


In  our  case.  rr  for  unstable  growth  must  be  of  the  order  of 

St  r- 


e  j  *_ 


In  terms  of  our  dimensionless  factors,  the  growth  K.s  -  277  S  /  >n  s,  is 
found  to  be 


Thus,  when  this  growth  is  larger  (and  hence  <  £  >  — >  2.  as  implied  by 
equation  (3.26)  ).  the  penetration  (diffusion)  distance  of  the  disturbances  is  much 
greater  than  the  propagation  wavelength.  This  means  that  magnetic  field  disturbances 
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will  be  minimal,  which  is  the  limit  postlated  in  our  earlier  simple  analysis.  Details 
are  plotted  in  figure  (3.2)  for  various  values  of  Vc  . 
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3.3  RESULTS 


First,  let  us  state  the  typical  background  situation  for  the  analysis. 


Electron  Temp. 

T 

-  20.000 

K 

Magnetic  field 

B 

-  0.075 

T 

Current  Density 

J 

-  5  x  105  Amp/m2 

Ionization  potential 

Vi 

-  15.8 

Volt. 

Electric  conductivity 

c r 

-  7212. 

A2/m2V 

So.  the  relations  between  those  non-dimensional  and  real  quantities  will  be. 

-  0.1291  x  10~3  u> 


V,  -  1.9094 

^  -  0.75  /  > 

The  wavelength  of  the  wave  is  of  the  order  of  1  mm  to  0.1  m.  since  the 
dimension  of  the  typical  MPD  devices  are  of  the  same  order.  ^  will  be  73  for 
*7\  being  1  cm.  and  so  the  range  of  -f  is  from  7  to  750. 

Also,  the  value  of  X  is  in  the  range  of  1  to  3.  As  appeared  in  equation 
(3.13).  it  depends  on  the  conductivity  >  .  magnetic  field  B  and  the  ionization 


potential  Vj.  Examining  the  typical  substances  being  used  to  generate  the  plasma.  Vj 
does  not  exceed  20.  In  order  to  achieve  higher  value  of  .  either  higher 


conductivity  or  lower  magnetic  field  is  required.  But  for  the  case  of  the  self-field 
MPD  devices,  the  value  of  B  is  related  to  the  current  density.  Also,  because  the 

plasma  is  m  Coloumb  collision  dominated  regime,  the  conductivity  depends  upon  the 
electron  temperature  and  thus  can  be  related  to  the  current  density,  in  other  words, 

the  parameter  is  a  function  of  current  density. 

**,  -  f  (  «  J  >  ) 

Since  the  onset  phenomena  limit  the  allowable  current  density.  fl" . 
ranges  only  from  1  to  3  for  the  existing  operation  regime. 

So.  the  only  relevant  region  of  figure  (3.2)  is  for  greater  than  49 

and  ranging  from  1  m  3.  i.e  the  flat  part  of  the  graph  where  the  critical  Hall 
parameter  is  nearly  2.  We  find  that  damping  due  to  magnetic  field  variations  does 
happen,  but  it  just  falls  out  of  the  range  of  the  operating  regime  of  the  MPD 

devices  currently  being  used.  Details  of  the  venations  of  "n  v*.  .  lfnag  with 

venous  <  fi  »  is  plotted  in  figure  (3.3). 

Looking  at  equation  (3.22) 


fie.  ' 


*  > 

Y. 


.  '~r  -  •*' 


(  .1,  — >  0  ) 
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We  can  mb  quite  easily  that  approaches  zero  for  both  limning  cases  where 

•4»  approaches  zero  or  infinity  This  agrees  with  our  finding  m  chapter  :»o  where 
there  is  only  damping  part  :n  the  propagation  of  disiurbacnas  when  we  'eiurr.  10  the 
simple  case  described  in  chapter  II  For  the  values  of  >  in  the  range  from  . 
mm  to  0.1  m.  is  plotted  w  wavelength  in  figure  (3  4)  It  is  growing  with  the 

wavelength  in  that  region  Looking  beck  to  our  numencai  eiprswom  regarding  the 
real  quantum  we  find  that  the  range  of  the  oscillating  pari  of  the  propagation  can 
be  up  to  "V  Hz 

In  figure  (3  5)  .*  imaginary  a  plotted  v*  the  curreat  angle  9  with 

vanoua  values  of  >  _»?  in  the  case  of  the  wavelength  being  a.21  centimeter  It  n 
apparent  that  $  -  135  is  the  moat  unstable  orients  nor  for  all  *  values  in 
this  case-  As  the  wavelength  gets  larger,  the  worst  orientation  becomes  less  than  135 
and  315  (Equation 
orientation  angle  a 
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3.23)  Bui  the  deviation  is  small  and  c  -  135  for  the  worst 
suit  an  acceptable  assumption 


3.4  CONCLUSIONS 


From  Uk  figures,  we  found  that  <  <3  >  -  2  is  still  the  threshold  value 
for  stability  at  the  worst  orientations  for  the  wavelength  from  1  mm  to  1  cm. 

The  variations  of  the  magnetic  field  has  little  effect  on  the  overall  stability 
picture.  In  the  existing  operation  regime  of  the  typical  MPD  device,  the  magnetic 
fwld  venations  appear  to  the  plasma  as  a  short  wave  effect  It  penetrates  throughout 
the  whole  region  quite  easily,  and  is  viewed  as  constant  inside  the  plasma,  because 
the  'skin -effect"  a  no  longer  effective  against  the  penetration.  So  the  assumption  of 
a  constant  magnetic  field  throughout  the  whole  plasma  stated  in  previous  chapter  is 
valid  and  sound. 


Also,  the  result  indicates  that  the  type  of  instability  is  still  the  same  as 
that  in  chapter  two.  where  it  is  caused  by  the  imbalance  between  the  energy  input 
(Joule  heating)  to  the  plasma  and  the  collisional  losses  by  the  electrons.  The  rate  of 
growth  of  the  instability  is  less  in  the  long  wavelength  case,  but  the  reduction  is  not 
enough  u>  reach  stability.  For  the  typical  MPD  devices,  whose  dimensions  are  of  the 
order  of  cenurneten  or  less,  wavelength  damping  appears  to  be  insignificant.  In  the 
next  chapter,  we  will  examine  the  effects  of  ionization  non-equilibrium  and  heat 
conduction  on  the  stability  picture. 
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CHAPTER  IV  PROPAGATION  OF  DISTURBANCES  THROUGH  THE  PLASMA 
OF  IONIZATION  NON-EQUILIBRIUM  WITH  VARYING  ELECTRON 
TEMPERATURE  BUT  CONSTANT  MAGNETIC  FIELD 

4.1  INTRODUCTION 

In  the  last  two  chapters,  we  examined  the  effects  caused  by  the  electron 
density  variation  by  looking  at  the  partially  ionized  regime  of  the  plasma.  We  found 
that  instability  occurs  when  the  energy  transfer  between  the  Joule  heating  and 
electron  collisional  lasses  is  not  in  balance.  The  magnetic  field  variation  was  found 
to  help  stabilize  the  situation  but  it  is  insufficient  in  the  regime  of  interest 

Now,  we  will  turn  to  the  situation  of  ionizaton  non-equilibrium.  In  other 
words,  the  kinetics  of  the  electrons,  rather  than  the  electron  energy  equation  will 
determine  the  ionization  rate(  Appendix  I  ).  The  excess  electron-ion  pairs  will 
diffuse  away  and  recombine  as  neutrals  at  the  boundry. 

Also,  instability  of  the  type  mentioned  in  the  last  chapter  does  not  appear 
near  full  ionization,  since  electron  density  can  no  longer  vary.  Now,  by  allowing  the 
electron  temperature  to  vary,  instabilities  might  occur  at  or  near  the  area  where  full 
ionization  is  reached. 

Also,  the  effect  of  heat  conduction  in  the  plasma  might  play  a  major  role 
in  the  overall  stability  picture.  In  previous  chapters  we  found  that  the  imbalance  of 
energy  transfer  causes  instability.  But  if  heat  conduction  occurs  within  the  plasma. 


4.2  CALCULATIONS 


4.2.1  ASSUMPTIONS 


As  shown  in  last  chapter,  the  effect  of  the  variations  of  the  magnetic  field 
is  insignificant  in  the  regime  of  our  interest  Thus,  we  will  return  to  the  assumption 
that  the  magnetic  field  remains  constant  (No  space  and  time  variations) 

The  positive  ion  and  neutral  temperatures  are  still  taken  to  be  constant  As 
stated  in  previous  chapters,  the  heavy  species  is  considered  as  immobile  when 
compared  with  the  electron  due  to  their  enormous  Compared  to  the  electron 
temperature  variation,  the  other  species’  temperatures  exhibit  only  smaller  variations. 
Therefore,  it  is  only  necessary  to  consider  the  electrons  in  the  number  conservation, 
momentum  and  energy  balances. 

So.  our  assumptions  will  be, 

1.  The  magnetic  field  is  constant.  (No  space  and  time  variations) 

2.  The  positive  ion  temperature  and  neutral  temperature  are  the 
same  and  remain  constant 

3.  In  the  background  state,  the  net  gain  in  electron  density  by  the 
ionization  process  will  be  balanced  by  the  am  bipolar  diffusion 
to  the  wall. 

The  heavy  species  number  density  (ion  density  plus  neutral 
density)  is  constant 


4. 


4.12  GOVERNING  EQUATIONS 


The  excess  energy,  which  will  fuel  the  ionization  process  and  transfer  away 
by  means  of  heat  conduction,  is  the  difference  between  the  Joule  heating  and  the 
collisional  losses  by  the  electrons.  So  the  energy  equation  will  be, 

*  'jf’  *  ^ )\  +  y  •  [»€Vlc  (*Ji  *  %  k  ~r  )J 

,,  (4.1) 

-1-7  *  4r*  -  3-^7  HcteCTc  -  ig  ) 


where  Tc  -  T.  £  "T*  Pc  ^  *£*- 

0  Veo 


The  convective  term  can  be  expressed  as 


*7  v<  i 


C 


r*  »  C 


(  n'c  will  be  suted  in  equation  (4.8)  ) 


The  heat  conduction  term  is  expressed  by. 


where  the  thermal  conductivity*  is 


~la(  *  Z3L  ^  T<  So)  jj  (■  k  Tt)  [  i  ^ 

Cl'1  pfTZT  *’  ua.  i  +  a* 

(4.3) 

m  Is 7  ‘if 3  (Te«,0~TjVj 
/  -t  /3X 


The  electron  momentum  equation  is, 


7  ->  T  *  /3  *■  O  £r  (4.4) 


where 


P>  m  /(*.  6 


c»t* 


The  plasma  is  still  assumed  to  be  in  Couloumb  collision  dominated  regime, 
so  the  conductivity  is  given  by  the  Spitzer-Hann  formula. 


7 


o  at i 

/  r\  /\ 


(4.5) 


,*  Sutton  ft  Sheman.  Engineering  Megnetohyd rodynemics,  pp  202 
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A.  - 

P„ 


X>o  »  J  ,  x  ,0 


Became  UA.  appears  as  a  weak  function  of  electron  temperature  and 

« 

density,  we.  as  in  previous  chapters,  assume  the  value  of  Ln  A  -  6. 

Maxwell’s  equations  are: 


"V  <  & 


“7  •  J  *  » 


Them  two  equations  lead  to  the  conclusion  that  Ey  and  Jx  are  constants,  in 
the  orientation  stated  in  previous  chapters. 


The  rate  of  change  of  electron  density  is  governed  by  the  kinetics  of 
electron  and  appeared  as. 


Oo 
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*CV,  '  V'VoV 


where 


n« 


_  =  5  =  e. 

Oo  v  K* 

is  the  Saha's  equation.  And  according  to  the  Bates'  law. 


Q ;  =  /  £>  9  *  i  & 


>  % 


The  convective  term  of  electron  density  variation  is  the  am  bipolar  diffusion 
to  the  wall. 


Da 

R* 


*1. 


(4.9) 


where  R  is  the  "effective  distance”  regarding  the  ambipolar  diffusion  process. 

To  relate  R  to  real  dimension,  let  us  model  the  plasma  in  the  following 
way:  It  is  confined  in  a  one  dimensional  space  with  length  H.  At  the  middle  of  the 
chamber,  there  is  no  electron  density  gradient  by  reason  of  symmetry.  Also,  electron 
density  at  the  wall  (  H/2  ;  -H/2  )  is  zero. 

Thus,  the  boundry  conditions  are: 


A 


*  o  ij-  % 


d  n, 

&<-• 

•j 


<=> 


r 
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The  diffusion  equation  is: 


2-  /  \  _  * 

d»  ^  sL>*  u«f  )  ~  (4.10) 

w/ 


Apply  the  boundry  conditions  and  we  find. 


»t>«. 


^  u*. 


(4.11) 


The  only  region  of  interest  for  the  analysis  of  this  chapter  is  the  middle 
region  of  the  model,  since  electron  density  is  diffusing  away  in  all  directions.  So. 


ncO}  *  H*  -- 

$  u>e>. 

and 

H  -  jT  R  (4.12) 

To  find  the  expression  for  the  diffusivity,  we  consider  the  am  bipolar 
momentum  balance  in  a  slab  geometry  similar  to  that  used  for  the  beat  diffusion 
calculation: 


1/+. 


J<(  ie  -»•  l«j  ) 


( r\z  +n~  , 


e 


(4.13) 
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and  so 


I 


» 

) 


rv 


^  ^  ^T<g  f  (a) 

L  P  5...  d'  J  'J 


Da.  the  diffusivity,  is  the  quantity  in  the  square  bracket. 


where 


^  1  *  xIUJ. 

f  c“ 


CL  frTV  ~ 
J  T/ 


(4.14) 


(4.15) 


/kin  is  the  reduced  mass  for  the  positive  ion-neutral  system.  Cjn  is  the 
effective  collisional  speed  of  positive  ions  with  neutrals,  and  Qjn  is  the  collisional 
cross-section  between  the  positive  ions  and  neutrals.  For  Argon.  Qjn  can  be 
represented  by  the  empirical  formula1 

Qin  -  19953  x  10’18  Tg“0-2 
(Qin  is  in  SI  unit  and  Tg  is  in  K) 


1  W.R  Cramer,  Journal  of  Chemical  Physics,  Vol.  30.  pp  641-642 
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Note  that  the  ambipolar  diffusion  rate  is  determined  by  the  positive 
ion-neutral  collisions  rate,  despite  the  fact  that  we  only  consider  the  electron  balance 


in  mass  conservation,  momentum  and  energy. 


The  electron  conservation  equation  is  then. 


k-  ^ 


'R.i  L 


-  'V  ) 


(4.16) 


The  electron  production  rate,  which  depends  on  electron  temperature  and 
on  the  departure  from  Saha  conditions  (as  the  plasma  is  assumed  in  ionization 
non-equilibrium),  will  equal  the  time  variation  of  the  electron  density  plus  the 
electron  flux  diffused  away  by  means  of  ambipolar  diffusion. 


4.13  CALCULATION  PROCESS 


The  coordinate  system  is  still  defined  as  before,  where  angle  5  is  the 
angle  between  the  wave  propagation  vector  and  the  average  current  vector.  (Fig.  2.1) 
Also,  a  plane  wave  dependence  of  the  disturbances  is  assumed. 


<u  > 


=  u  e 


The  disturbances  will  exhibit  unstable  behavior  if  the  imaginary  part  of 
u>  is  negative  (  M  is  real) 


Now,  let  us  define  some  non-dimensional  parameters  befcre  we  arrange  the 
governing  equations. 


(4.17) 


is  the  ratio  between  the  thermal  energy  of  the  electrons  to  the  energy  required  for 
the  ionization  process. 


(4.18) 


is  the  rate  of  electron  loss  by  means  of  ambipolar  diffusion.  (It  has  a  dimensional 


(4  19) 


-  ~  0  i  - 


■  y 


And  finally. 


d  ~ 


n* 


(4.20) 


the  degree  of  ionization.  (  npj  -  n«  +  nn  is  the  heavy  species  density  ) 

Applying  the  perturbation  treatment  (same  as  those  in  chapter  three)  to  the 
separate  terms  of  the  energy  equation  (4.1),  we  find  the  following  terms  to  be: 

Time  variation  term: 


z  t 

v  T  K  Tc  }  'e  |  i  t  'S  <av>  (  -*  f  l«TTc.  >  )  ? 


cO  <  / 1  ? ,  ~  '  "  r  >  ofe 


Convective  term: 


^7  •  L  >1*Vc  *  t 


^  cl/ 1  (  /  !  C-d  ~  Z  ) 

+  >*'**''■ ■'] 

+  -v  <  >>  49  -+  <  ft  -r  ^<7,;  '♦U 

"  (  ''A  (itirJ  [  t'  ~  <£.  •></’ftv  <- ^ 

-  ai  <',<,><’»  lC  >  J 

—  t  >  pvA  s  -  cr  .  | 
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Heat  conduction  term: 


*7  ■  A*'  — 


,  i 

13  7  .  ^  3  (^'Vio fij 


I  -  «i> 


Joule  heating  term: 


<7* 


=  -  ^-ColcIS^  V  '•  ^  \3?  Sfidb  <5/i* 


Collisional  loss  term: 


-e.  - 


-  —  <  t*  >  «■-'* 


eA  ^  St  *>  S>  ’c 


And  the  energy  equation  becomes: 


»V  ^Vt  (/-•  <jr  j  ('■*  -  i  j  J£»  ■«  >  >.  ■  r  >  J 


*4.  'c  >-  -+  iH'k'-/  'rss- 


/i  7  V Li  <Vc’ 


v* 

y  <  Tr  > 


'I 


■+  — 


V>-  -  - 


7  <tf> 


I 


/  _  > 

A 


(4.21) 


r* 

{^vf 


-  Jhf  -  'X.><vc>  v->  j 


Y~7  </>-•—&  Smc.  ~ 
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Reamnfing  the  mas  conservation  equation,  we  gee 


Convective  term: 

c  kr  (  i  ) 


Time  dependent  term: 


^•L 


-  <‘0  </!.<•■>  Sk^ 


Total  growth  rate: 


l«.  -  ‘  z  )  *  *?  'V  ><»ie>  j 

-t  ^  A',V  'J  >  v '  >”  >  J  <5,tV 


Overall  conservation  equation: 


s  „  *  /  -  » 

-♦  <*•  -*<w>  •  e>  v  el't 


{  ^  (  at  W 


< ,  -><v>e  V 


/  *;  -  * 


(4.22) 
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Putting  the  governing  equations  into  matrix  form,  we  get. 


) 


<3 


where 


> 

I 


A(  1.  1  )  -  i  e  J,<0. e'tz  *  ?\/t  (l-rZt  >[ 

-r  *f<?.'><W>  ]  •+  Te  <& 

.  £  ✓.  > 

-»  i  /<  us  — 

e 

-»  <7~f>  ^  * 

~f  -r  4^  ►  -  3l  v  ft  > 


A(  1.  2  )  - 


I 

I 


« i/.  c  i  ~  2, .[  -±  -  <\\ > <  ?  o' 


A(  2,i  } 


■fi  (  ^  - «?  ; 

—  A  <  >  fV{*>  ■2* 


I  A(  2,  2  )  - 


(OOe>  -*  -Jt'  <?,'  ■»  <  Opy 


(4.22  b) 


(4.22  c) 


(4.22  d) 


(4.22  e) 


(4.22  f) 


cO.  the  variable  of  interest,  can  be  solved  for  by  setting  the  determinant  to 


zero. 


Other  parameters  such  as  the  electron  density,  current,  etc.,  are  calculated 
as  follow:  From  the  zeroth  order  continuity  equation. 


<■*; ~ (4.23) 


Rearranging  it  and  we  find. 


<*<>  _ 


_ _ <  -i 


(4.24) 


where  njj  and  are  input  parameters.  The  rest  just  shows  electron  temperature 
dependence.  Figure  (4.1)  shows  the  relationship  between  electron  density  and 
temperature  for  a  particular  density.  Other  parameters  such  as  pressure  and  input 
current  will  be  shown  in  Table  (4.1). 

Then,  the  energy  loss  by  means  of  collisions  between  electrons  and  ions 
can  be  calculated  as. 


<&> 


k.<  Te> 


(4.25) 


From  the  zeroth  order  energy  equation,  we  find  the  input  current  to  be. 


(4.26) 


Then.  can  be  calculated  from  the  full  matrix,  and  will  be  presented  in 


the  next  section. 
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4.3  RESULTS 


Solving  the  full  determinant  equation  (4.22  b).  we  find  in  general  two 


roots  for  (  u>r  +  *  ui  j  )  for  each  set  of  M .  fi  .  o(  ,  etc.  This  means  that  there 


are  two  distinct  modes  of  plasma  fluctuation  (and  possibly  growth).  In  order  to 


understant  these  two  modes,  we  first  examine  the  limiting  conditions  where  one  or 


the  other  will  predominate. 


First,  very  near  full  ionization  (for  *  »  0.99.  as  we  will  see),  electron 


density  variation  becomes  negligible,  and  the  energy  equation  controls  the  growth. 


With  reference  to  equation  (4.22  c).  we  see  that  9  -  90*  is  the  worst  orientation. 


and  that,  since  at  9  -  90*  all  terms  in  equation  (4.22  c)  are  real.  «*>  will  turn  out 


purely  imaginary  (i.e..  the  growth  or  decay  will  be  purely  local).  Also  the  Saha 


factor  <  S  >  will  be  very  large,  implying  very  fast  electron  kinetics,  so  that  Saha 


equilibrium  will  be  maintained  and  the  perturbation  continuity  equation  becomes. 


<«?,>  W'T  A 


-O?-  ><  s>  <-n«> 


(4.27) 
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The  energy  equation  becomes. 


<'rv>'V«  e'j.  ■+  H*s<.  Co*^u>-- 
—  Z  <Br>  •*  */■  >§  <fl.  *  *',*»><  >"j 


"  (■  ;  y  2*  )  <  /?**>  «'/7e ->*"«■  -  ■>  &* 


_ 


*■  c> 


(4.28) 


Combining  the  two  expressions  above  we  get. 


—  ~  *■  S'  ">  "*  <A7 'v  S  >  «  >> 


-  p/-  f  -  <  ^  N  < >\»>  <■  s >  .  ■:  j<_-‘_^.-. 


(4.29) 


and  so.  « ^  ?c  <  ,u>  Ot  ~  A/*  r-  ;  <r.  ■»•.  -  f  r-*^ 

—  j*  <  t*  > 

The  plasma  is  unstable  without  heat  conduction  and  this  is  a  static 
instability,  in  the  sew  that  no  wave  propagation  is  involved.  To  look  for  the  Hall 
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parameter  under  which  condition  the  plasma  with  heet  conduction  is  at  the  neutral 
stability  point,  we  set  ^  equal  to  zero  for  the  above  equation. 


H* 


Cent  «.  f 


,Zl  Vjx.o*  c  <7’>  (*>*)}' 

~*\1’  (  i  *  <  3  /  > 


(4.30) 


*  fir  > 


More  accurately,  we  set  <9  -  90  *  and  u>.  •  0  for  the  full  matrix  to 
solve  for  the  neutral  Hall  parameter.  This  is  plotted  in  figure  (4.2)  for  various 
values  of  electron  temperature.  It  shows  that  there  is  a  stabilizing  effect  of  beat 
conduction  within  the  plasma,  since  at  wavelengths  greater  than  about  3  cm.  the 
plasma  will  be  unstable  at  tI  >  99%.  whereas  it  will  not  for  shorter  waves. 

In  the  opposite  limit,  where  electron  density  variation  is  the  dominant 
effect  (low  degree  of  ionization),  we  combine  the  two  governing  equations  with 
Six  "  o. 


t  ,-><-«>  ><->>  1-  <  t\  ><\S  S' -  O 


I 

t 

1  and 


jjjj  r  t 

x  ^  >  er  e/  It-  l,)  ^  <  ?.  s  KS  >  <  •’•»> 

■*  Jx  -  if  ( 
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*  : 


Since  in  this  limit  electron  kinetics  must  be  irrelevant,  we  eliminate  <  Ri 


> 


between  the  two  equations.  The  damping  pert  of  ^  is  then 


U\ 


( 


>'v 


>5 


<£. 


(4.31) 


The  worst  orientation  returns  to  &  -  135 and  the  stability  factor  is 
once  again  the  Hall  parameter.  The  plasma  behavior  returns  to  the  same  picture 
diacuaed  in  the  previous  chapters,  except  that  the  critical  Hall  parameter  is  smaller 
in  order  to  account  for  the  increased  Joule  dissipation  due  to  the  ambipolar  diffusion 
effect  (Note  that  in  the  regime  where  the  plasma  is  near  full  ionization,  this  limit 
will  not  be  observed,  and  the  critical  Hall  parameter  has  to  be  found  by  solving  the 
full  matrix.) 

By  solving  the  determinant  above,  we  achieve  two  solutions  of  the  damping 
part  of  oo  ,  which  correspond  to  the  two  modes  of  stability  described  above. 
Then,  by  concentrating  on  the  worst  angles,  we  can  map  the  behavior  of  as 

the  degree  of  ionization  changes.  As  expected,  the  electron  temperature  variation 
mode  shows  up  near  the  full  ionization  regime  (Figure  (4.3)),  while  the  density 
variation  mode  dominates  the  rest  (Figure  (4.4)). 


Stressing  our  interest  in  the  situation  where  the  plasma  is  in  neutral 
stability  for  the  worst  orientation,  we  examine  the  ratio  between  *te  and  Snt  : 

For  the  solution  corresponding  to  the  temperature  variation  mode,  'sCt  is 
the  dominant  factor  for  stability  analysis.  It  only  appears  in  the  region  where  the 
degree  of  ionization  is  very  close  to  unity,  and  it  agrees  with  the  description  stated 
previously.  The  neutral  Hall  parameter  vs.  degree  of  ionization  is  plotted  in  figure 
(4.5)  (Same  as  the  X  -  0.02  m  curve  in  figure  (4.2)  ). 

For  the  other  solution,  corresponding  to  the  density  variation  mode,  Stic  is 
the  dominant*  factor  where  the  plasma  is  partially  ionized  and  its  contribution 
weakens  as  the  degree  of  ionization  rises  (Figure  4.6).  At  the  regime  where  both 
variations  arc  comparable  as  one  mode  begins  to  diminish  and  the  other  to  rise,  the 
coupling  effects  between  the  two  variation  modes  become  important;  the  critical  Hall 
parameter  appears  to  be  the  lowest  (Figure  4.7),  and  wave  propagation  (as  measured 
by  the  real  part  of  ui  )  is  diminishing  as  the  frequency  approaches  zero  (Figure 
4.8).  As  the  plasma  is  further  ionized,  Sne  variation  becomes  insignificant  The 
critical  Hall  parameter  rises  again,  and  wave  propagation  resumes. 

The  overall  boundary  of  the  Hall  parameter  for  stability  regardless  of  the 
mode  of  variation  is  plotted  in  figure  (4.9),  for  a  2  cm  wavelength.  It  shows  a 
smooth  transition  between  the  two  modes. 


4.4  CONCLUSION 


We  have  explored  the  effects  of  heat  conduction,  which  indeed  is  a 
stabilizing  agent  to  the  plasma.  Since  the  direction  for  the  conduction  is  assumed  to 
be  transversal  to  the  magnetic  field,  the  thermal  conductivity  depends  on  the  Hall 
parameter  inversely.  So  at  high  Hall  parameter,  this  effect  diminishes  and  can  no 
longer  provide  the  stabilizing  force.  The  associated  wavelength,  must  be  in  the  order 
of  0.01  meter  for  this  stabilization  to  be  visible.  Also,  this  mechanism  affects  the 
temperature  variation  mode  only.  For  the  other  mode,  the  models  from  previous 
chapters  remains  the  proper  descriptions. 

Electron  temperature  variation  leads  to  the  regime  of  static  instability.  It 
only  happens  at  the  regime  where  density  variation  can  no  longer  be  achieved,  i.e. 
the  plasma  is  almost  fully  ionized.  At  lower  ionization  fractions,  this  instability  does 
not  appear,  and  the  density  variation  becomes  the  dominant  factor  for  stability 
concern.  The  coupling  effect  between  the  two  modes  makes  the  transition  between 
the  regimes  smoothly,  but  the  stability  of  the  £n«  mode  deteriorates  markedly  in  the 
region  where  both  modes  interact,  i.e.,  between  of  0.9  and  0.99.  This  is  a  regime 
of  interest  for  MPD  thrusters,  and  more  study  should  be  devoted  to  it. 
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CHAPTER  V  SUMMARY  AND  CONCLUSIONS 


We  have,  throughout  the  course  of  this  thesis,  established  some  microscopic 
physical  models  for  the  plasma  contained  in  a  co-axial  chamber  with  current  passing 
through  between  the  two  end  walls,  so  we  can  examine  various  effects,  such  as 
electron  temperature,  density  and  magnetic  field  variations,  etc.,  that  affect  the 
overall  local  stability  picture.  The  following  conclusions  have  been  drawn: 

For  a  partially  ionized  plasma  of  constant  induced  magnetic  field, 
instabilities  occur  whenever  the  direction  of  the  wave  vector  of  the  propagating 
disturbances  falls  outside  the  direction  of  the  current  flow  and  electric  field,  since 
Joule  heating  in  the  originally  low  density  region  will  decrease  further,  and  thus  it 
will  increase  the  discrepancy  of  electron  density  between  different  regions.  So  the 
uniform  background  can  no  longer  be  attained. 

The  threshold  value  for  the  Hall  parameter  for  which  the  disturbances  are 
in  neutral  stability,  i.e.  no  growth  or  damping  (critical  Hall  parameter),  is  2  in  the 
worst  wave  vector  orientation. 

Even  though  the  induced  magnetic  field  exhabits  space  and  time  dependence 
associated  with  the  presence  of  currents,  the  above  threshold  value  of  the  Hall 
parameter  is  still  valid  for  the  wavelengths  of  spatial  variation  of  the  magnetic  field 
from  1  mm  to  1  cm.  This  variation  has  only  little  effects  on  the  overall  stability 
picture.  The  spatial  variation  of  the  magnetic  field  does  provide  some  stablizing 
effects,  where  the  growth  rate  of  instabilities  is  less  for  a  long  wavelength.  But  this 
damping  is  insufficient  for  the  typical  MPD  devices.  The  magnetic  field  variations 


penetrate  throughout  the  whole  region  quickly  .  And  so  it  is  viewed  as  constant 

inside  the  plasma,  while  the  "skin-effect"  which  acts  against  the  penetration  is  not 
effective  enough.  Constant  magnetic  field  throughout  the  plasma  is  a  valid 
assumption  for  those  operations. 

For  a  fixed  magnetic  field,  two  potentially  unstable  modes  of  wave 
propagation  or  growth  exist.  These  are  characterized  by  the  predominance  of  electron 
density  or  electron  temperature  fluctuation  in  each  of  them.  The  mode  with 
predominant  density  fluctuation  tends  to  be  the  main  instability  at  low  ionization 
fraction,  while  the  electron  temperature  mode  becomes  unstable  near  full  ionization 
(  d  ■  0.99).  Of  course,  both  modes  occur  throughout  the  range  of  ionization 
fraction.  No  wave  propagation  is  invloved  in  this  electron  temperature  mode 

instability.  As  the  ionization  fraction  decreases,  it  diminishes  very  rapidly. 

Heat  conduction  is  a  stabilizing  agent  to  the  undesirable  situation  described 
above.  But  as  the  direction  for  the  heat  conduction  is  transversal  to  the  magnetic 
field,  the  corresponding  thermal  conductivity  depends  on  the  Hall  parameter,  and  so 
this  effect  is  insufficient  to  alter  the  stability  picture  at  high  Hall  parameters.  The 
corresponding  wavelength  for  the  heat  conduction  must  be  of  the  order  of  0.01 
meter  or  less  in  order  for  this  stabilizing  effect  to  be  effective. 

Because  of  the  coupling  effect,  a  smooth  transition  can  occurs  between  the 

electron  density  and  temperature  variations  in  the  region  (  ^  -  0.9  to  0.99)  where 

both  modes  interact  with  each  other.  But  this  effect  also  courses  the  stability  of  the 
electron  density  variation  mode  to  deteriorate  to  an  undesirable  level. 


Finally,  since  the  operating  regime  of  the  typical  MPD  thrusters  falls  into 
the  region  of  the  coupling  between  the  two  variations  for  the  electrons,  it  will  be 
necessary  to  devote  more  study  into  the  coupling  effect  between  the  two  modes. 
Other  effect  such  as  variations  of  plasma  properties  along  the  propagating  wave  front 
might  affect  the  overall  stability  picture,  and  it  is  also  worthwhile  to  investigate. 
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Figure  1.1 

Seif-field  coaxial  type  accelerator 


Figure  1.2 

Self-field  coaxial  type  accelerator 
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Figure  11 

Coordinate  system  for  the  analysis 
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Figure  14 

Diagram  of  averaged  current  direction  in  stable  regime 


Figure  15 

Diagram  of  averaged  current  direction  in  unstable  regime 


Figure  3.1 

Variation  of  the  wont  current  angle  vs.  wavelength 
ranging  from  6.28  mm  to  0.628  m 


Figure  3.2 

Critic*!  Hall  parameter  vs.  wave  number 
for  spatial  magnetic  field  variation 
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Degree  of  Ionization 


Figure  4.1 

Degree  of  ionizttion  vs.  electron  temperature 


Density  :  0.3345  x  10‘Vg/M3 
Plasma  Dimension  *  0.05  M 


Neutral  Hal  I  Parameter 


Figure  4.2 

Neutral  Hall  parameter  vs.  degree  of  ionization 
for  various  values  of  wavelength 
(Sl2  mode) 


Density  1  0.3345  x  10  3 Kg/M2 
Plasma  Dimension  t  0.05M 
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Figure  4.3 


Damping  pert  of  U>  vs.  degree  of  ionization 
(^4  mode) 


Hal  I  Parameter  :  3.5 
Density  *  0.3345  x  10"3  Kg/M3 
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Figure  4.4 

Damping  part  of  ^  vs.  degree  of  ionization 

mode) 

Hall  Parameter  *•  3.5 
Density  :  0.3345  x  10  1  Kg/M3 
Plasma  Dimension  ’  0.05  M 
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Figure  4.5 

Neutral  Hall  perimeter  vs.  degree  of  ionization 
(£4  mode) 


Wavelength  :  0.2  x  10'1  M 
Density  s  0.3345  x  lO'^Kg/M 
Plasma  Dimension  t  0.05  M 
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Wave  Frequency 


Figure  4.8 

Wave  frequency  vs.  degree  of  ionization 
(§n^  mode) 


Wavelength  *  0.2  x  10''  M 
Density  *  0.3345  x  10*3  Kg/M3 


Degree  of  Ionization 
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«#.«•  e.as.-i  1.25*0 


Figure  4.9 

Hall  parameter  vs.  degree  of  ionization 
(hn«  and  modes) 


Wavelength  *  0.2  x  10 '  M 
Density  ••  0.3345  x  10 3  Kg/M3 
Plasma  Dimension  *•  0.05  M 
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Table  4.1 


Pressure  and  input  current  for  various 
degrees  of  ionization 


Gas  temperature 

Heavy  species  number  density 
Plasma  dimension 

:  1000 ’K 
:  5  x  10 
:  0.05  Meter 

Electron 

temoerature 

(xlO5  K) 

Degree  of 
ionization 

Pressure 

(N  /  M') 

Current 

(x  107  J/Ma  ) 

.1000 

.1928 

202.1 

0.39786-01 

.1100 

•  1.394 

1*02.7 

0.9168E-01 

.1150 

•  5758 

526.2 

.1221 

.1200 

•  7006 

61*9.1* 

.1513 

.1250 

.8009 

760.1 

.1765 

.1300 

.8728 

852.3 

.1963 

.1350 

.9202 

926.6 

.2112 

.uoo 

•  91*99 

987-1 

.2225 

.1500 

.9793 

1083- 

.2386 

.1600 

•  9906 

1163- 

.2505 

.1800 

•  9975 

1309. 

.2702 
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APPENDIX  I 


SAHA’S  EQUATION  AND  EQUATION  FOR  IONIZATION  AND 
RECOMBINATION  RATE 

For  an  ionization  process, 

t  z-  -*■  <■  ^  O  (app  1.1) 

where  z  is  the  number  of  electron  released  if  one  positive  ion  is  produced.  e.i.n 
represent  the  electron,  ion  and  neutral  respectively. 

From  the  law  of  Mass  Action,  we  can  write 

( Vg.;  (Ay,  )  (&>e  2 

■/.  )  frPP 

Ne  :  number  of  electron 

Qe  :  partition  function  for  electron 

Nj.  Nq  Qj  .  Qn  are  corresponding  functions  for  ion  and  neutral(atom). 

(  z  -  1  for  Agron  plasma,  which  is  assumed  to  be  used  in  our  study.  ) 


The  partition  functions  for  ail  three  species  are. 


/? 


kt 


oV-  0  for  electron  and  ion  since  they  don’t  possess  any  energy  at  rest 

£  -  -eVj  for  neutrals  at  rest.  (It  is  also  called  the  energy  for  ionization) 
Vj  is  the  ionization  potential(voltage). 


+r 

€ 


:/ 


SiC 


s>. 


-A 

y  v  (^5  /  *»-  c. 


) 


^  ^  T’a 


o 


Qnu»t  is  the  internal  partition  function. 


Then,  equation  (app  1.1)  becomes. 


determined  only  by  temperature. 


Now,  if  a  net  ionization  rate  does  exist,  then  we  can  write. 


(  n^  )  -  (  n^  ionization  “  (  &e  ^recom  bination 


Assuming  ionization  is  by  electron-neutral  impact,  and  recombination  is  by 
a  three-body  e-e-i  process. 


(  vn'e  )  =  S'  rT«  -  o(  nt  o;  <*PP  1.5) 


(  if  .  d.  are  function  of  electron  temperature  only.  ) 

In  equilibrium,  nc  •  0  (Let  us  denote  r>e  as  the  electron  density  at 
equilibrium) 


T\  *  J  4  \ 

where  ^  ,  is  the  Saha  equilibrium  constant  depends  on  temperature  only. 


Return  to  (app  l.S),  we  find 


^e. 


ck 


(app  1.6) 
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where  £  is  expressed  by  the  Bates'  law. 


Jt  -  1.09  x  10-20  Te  -4-5 

If  the  Saha  equilibrium  density  (  )  is  much  greater  than 

the  electron  density(plasma  density),  then 


r>^ 


in  other  words,  the  bulk  recombination  process  is  negligible.  In  a  steady  state  this  is 
usually  due  to  am  bipolar  diffusion  to  walls,  followed  by  recombination  there.  In  a 
transient  state,  it  may  occur  at  the  start  of  a  strong  discharge,  before  the  electron 
density  has  had  time  to  build  up. 
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